Wwff 

non 


US  Army  Corps 
of  Engineers® 

BVP^^Sli?8tory 


Two-Dimensional  Analysis  of  Natural 
Convection  and  Radiation  in  Utilidors 


Paul  W.  Richmond 


June  1999 


IHM 


^proved  for  PuilpMease 
tetlsibufiotfflBli" 


19990720  003 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestion  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington, 
VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 


1 .  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE 

_  June  1999 


4.  TITLE  AND  SUBTITLE 

Two-Dimensional  Analysis  of  Natural  Convection  and 
Radiation  in  Utilidors 


6.  AUTHORS 

Paul  W.  Richmond 


3.  REPORT  TYPE  AND  DATES  COVERED 


5.  FUNDING  NUMBERS 

PE:  4A762784AT42 
PR:  6.27.84A 
TA:  M 
WU:  07 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory 
72  Lyme  Road 

Hanover,  New  Hampshire  03755-1290 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

CRREL  Report  99-7 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  the  Chief  of  Engineers  DoD  HPC  Center 

Washington,  DC  20314-1000  Arctic  Region  Supercomputing  Center 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


1 1 .  SUPPLEMENTARY  NOTES 

For  conversion  of  SI  units  to  non-SI  units  of  measurement,  consult  ASTM  Standard  E380-93,  Standard  Practice  for  Use  of  the 
International  System  of  Units,  published  by  the  American  Society  for  Testing  and  Materials,  1916  Race  St.,  Philadelphia,  Pa. 
19103. 


12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  is  unlimited. 
Available  from  NTIS,  Springfield,  Virginia  22161 


A13.  ABSTRACT  {Maximum  200  words) 

Central  heating  plants  are  often  used  on  large  building  complexes  such  as  university  campuses  or  military  bases. 
Utilidors  can  be  used  to  contain  heat  distribution  lines  and  other  utilities  between  a  utility  station  and  serviced 
buildings.  Traditional  thermal  analysis  of  utilidors  is  one-dimensional,  with  heat  transfer  correlations  used  to 
estimate  the  effects  of  convection,  radiation,  and  two-dimensional  geometric  effects.  The  expanding  capabilities 
of  computers  and  numerical  methods  suggest  that  more  detailed  analysis  and  possibly  more  energy-efficient 
designs  could  be  obtained.  This  work  examines  current  methods  of  estimating  the  convection  and  radiation  that 
occur  across  an  air  space  in  square  and  rectangular  enclosures  and  compares  them  with  numerical  and  experi¬ 
mental  data. 


14.  subject  terms  Central  heat  distribution 
Effective  conductivity 
Enclosures 

Numerical  heat  transfer 


Radiation  and  convection 

Utilidors 

Utilities 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 


UNCLASSIFIED 


NSN  7540-01-280-5500  “L 


UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


one  QUALITY  INSPECTED  4 


15.  NUMBER  OF  PAGES 

78 


16.  PRICE  CODE 


20.  LIMITATION  OF  ABSTRACT 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  Z39-18 
298-102 


Abstract:  Central  heating  plants  are  often  used  on  large 
building  complexes  such  as  university  campuses  or 
military  bases.  Utilidors  can  be  used  to  contain  heat 
distribution  lines  and  other  utilities  between  a  utility  sta¬ 
tion  and  serviced  buildings.  Traditional  thermal  analy¬ 
sis  of  utilidors  is  one-dimensional,  with  heat  transfer 
correlations  used  to  estimate  the  effects  of  convection, 
radiation,  and  two-dimensional  geometric  effects.  The 


expanding  capabilities  of  computers  and  numerical 
methods  suggest  that  more  detailed  analysis  and  pos¬ 
sibly  more  energy-efficient  designs  could  be  obtained. 
This  work  examines  current  methods  of  estimating  the 
convection  and  radiation  that  occur  across  an  air  space 
in  square  and  rectangular  enclosures  and  compares 
them  with  numerical  and  experimental  data. 


Cover:  Construction  of  a  utilidor  at  Eielson  AFB,  Alaska.  (Photo  courtesy  of  Dr.  Gary  Phetteplace.) 


How  to  get  copies  of  CRREL  technical  publications: 

Department  of  Defense  personnel  and  contractors  may  order  reports  through  the  Defense  Technical  Information  Center: 
DTIC-BR  SUITE  0944 
8725  JOHN  J  KINGMAN  RD 
FT  BELVOIRVA  22060-621 8 
Telephone  1  800  225  3842 
E-mail  help@dtic.mil 

msorders@dtic.mil 
WWW  http://www.dtic.mil/ 

All  others  may  order  reports  through  the  National  Technical  Information  Service: 

NTIS 

5285  PORT  ROYAL  RD 
SPRINGFIELD  VA  22161 

Telephone  1  800  553  6847  or  1  703  605  6000 

1  703  487  4639  (TDD  for  the  hearing-impaired) 

E-mail  orders@ntis.fedworld.gov 

WWW  http://www.ntis.ciov 

A  complete  list  of  all  CRREL  technical  publications  is  available  from 

USACRREL  (CEERD-IM-HL) 

72  LYME  RD 

HANOVER  NH  03755-1290 

Telephone  1  603  646  4338 

E-mail  techpubs@crrel.usace.army.mil 

For  information  on  all  aspects  of  the  Cold  Regions  Research  and  Engineering  Laboratory,  visit  our  World  Wide  Web  site: 
http://www.crrel.usace.army.mil _ 


Go  to  Contents 


CRREL  Report  99-7 


Two-Dimensional  Analysis  of  Natural 
Convection  and  Radiation  in  Utilidors 

Paul  W.  Richmond 


Prepared  for 

OFFICE  OF  THE  CHIEF  OF  ENGINEERS 

Approved  for  public  release;  distribution  Is  unlimited. 

Back  to  Contents 


June  1999 


PREFACE 


This  report  was  prepared  by  Dr.  Paul  W.  Richmond,  Mechanical  Engineer, 
Applied  Research  Division,  U.S.  Army  Cold  Regions  Research  and  Engineering 
Laboratory  (CRREL),  Hanover,  New  Hampshire.  Funding  for  this  effort  was  pro¬ 
vided  by  DA  Project  4A762784AT42,  Cold  Regions  Engineering  Technology ,  Work  Unit 
M07,  Heat  Transfer  from  Buried  Utility  Lines .  The  work  was  also  supported  in  part 
by  a  grant  of  HPC  time  from  the  Department  of  Defense  HPC  Center,  Arctic 
Region  Supercomputing  Center  (denali.edu). 

The  author  expresses  his  appreciation  to  the  people  at  CRREL  who  assisted  in 
this  project.  Special  thanks  go  to  Ronald  Farr,  who  constructed  the  experimental 
apparatuses.  The  author  particularly  thanks  Dr.  Ronald  Liston  for  his  continued 
encouragement  and  support.  This  report  was  reviewed  by  Dr.  John  P.  Zarling,  Dr. 
Thomas  Kinney,  Dr.  Debendra  K.  Das,  and  Dr.  Gary  Gislason. 

The  contents  of  this  report  are  not  to  be  used  for  advertising  or  promotional 
purposes.  Citation  of  brand  names  does  not  constitute  an  official  endorsement  or 
approval  of  the  use  of  such  commercial  products. 


ii 


Go  to  Contents 


CONTENTS 


Page 

Preface .  u 

Nomenclature  and  abbreviations .  vi 

Introduction .  1 

Background  ... . 3 

Numerical  model .  13 

Interpolation  functions .  16 

Solution  procedure .  17 

Radiation .  20 

Model  verification .  21 

Experimental  procedure .  24 

Experimental  apparatuses . 24 

Data  acquisition  system .  26 

Procedure .  29 

Results .  30 

General  information . 30 

Experimental  data .  32 

Numerical  data .  35 

Analysis .  39 

Summary . 49 

Conclusions  and  recommendations .  49 

Literature  cited .  50 

Appendix  A:  Experimental  results . 53 

Appendix  B:  Numerical  results .  57 

Abstract .  68 


ILLUSTRATIONS 

Figure 

1.  Cross  sections  of  two  utilidors  constructed  in  the  Arctic .  2 

2.  Heat  transfer  correlations  for  vertical  enclosures .  5 

3.  Cross  section  of  a  concentric  pipe  conduit .  7 

4.  Rectangular  enclosure  configuration  of  Ghaddar .  8 

5.  Experimental  configuration  of  Stewart  and  Verhulst .  9 

6.  Experimental  apparatus  of  Babus'Haq  et  al .  10 

7.  Heat  transfer  equations  for  pipes  in  enclosures .  10 

iii 


Back  to  Contents 


Figure  Page 

8.  Generic  utilidors  for  current  utilidor  thermal  analysis  procedure .  11 

9.  Curved  isoparametric  quadrilateral  element .  16 

10.  Viewfactor  analysis  of  FlH 6 .  20 

11 .  Two-dimensional  conduction  problem .  22 

12.  FECOME  results  for  a  vertical  square  enclosure  at  Rayleigh 

number  105 .  23 

13.  Grid  used  for  comparing  viewfactors  and  radiation  heat  flux 

calculations .  23 

14.  1-ft  x  1-ft  experimental  apparatus  without  lid .  24 

15.  Thermal  conductivity  of  expanded  polystyrene .  25 

16.  Schematic  cross  section  of  the  1-ft  x  1-ft  apparatus .  25 

17.  Thermal  conductivity  of  fiberglass  pipe  insulation .  26 

18.  2-ft  x  4-ft  enclosure  with  4-in.  and  8-in.  insulated  pipes .  27 

19.  Schematic  diagram  of  the  2-ft  x  4-ft  enclosure .  27 

20.  Data  acquisition  system  for  use  with  the  1-ft  x  1-ft  enclosure .  29 

21.  Nusselt  and  Rayleigh  number  plot  for  the  4-in.  pipe  in  the  1-ft  x  1-ft 

enclosure  (experimental  data) .  33 

22.  Nusselt  and  Rayleigh  number  plot  for  the  insulated  4-in.  pipe  in 

the  1-ft  x  1-ft  enclosure  (experimental  data) .  33 

23.  Nusselt  and  Rayleigh  number  plot  for  the  2-in.  pipe  in  the 

1-ft  x  1-ft  enclosure  (experimental  data) .  34 

24.  Nusselt  and  Rayleigh  number  plot  for  the  2-in.  insulated 

pipe  in  the  1-ft  x  1-ft  enclosure  (experimental  data) .  34 

25.  Nusselt  and  Rayleigh  number  plot  for  the  2-ft  x  4-ft  enclosure 

(experimental  data) .  35 

26.  Mesh  for  the  1-ft  x  1-ft  enclosure  with  a  4-in.  pipe .  35 

27.  Mesh  for  the  2-ft  x  2-ft  enclosure  with  a  2-in.  insulated  pipe .  36 

28.  Mesh  for  the  2-ft  x  4-ft  enclosure,  4-in.  and  8-in.  insulated  pipes .  36 

29.  Mesh  for  the  2-ft  x  2-ft  enclosure  with  two  2-in.  pipes .  36 

30.  Nusselt  and  Rayleigh  number  plot  for  the  4-in.  pipe  in  the 

1-ft  x  1-ft  enclosure  (numerical  data) .  37 

31.  Nusselt  and  Rayleigh  number  plot  for  the  insulated  4-in.  pipe 

in  the  1-ft  x  1-ft  enclosure  (numerical  data) . .  37 

32.  Nusselt  and  Rayleigh  number  plot  for  the  2-in.  pipe  in  the 

1-ft  x  1-ft  enclosure  (numerical  data) .  38 

33.  Nusselt  and  Rayleigh  number  plot  for  the  insulated  2-in.  pipe 

in  the  1-ft  x  1-ft  enclosure  (numerical  data) .  38 

34.  Nusselt  and  Rayleigh  number  plot  for  the  1.27-ft  x  1.27-ft 

and  2-ft  x  2-ft  enclosures  (numerical  data) .  38 

35.  Heat  flux  from  the  4-in.  pipe  through  the  1-ft  x  1-ft  enclosure .  39 

36.  Heat  flux  from  the  2-in.  pipe  through  the  1-ft  x  1-ft  enclosure .  40 

37.  Heat  flux  from  the  4-in.  insulated  pipe  through  the  1-ft  x  1-ft 

enclosure .  40 

38.  Heat  flux  from  the  2-in.  insulated  pipe  through  the  1-ft  x  1-ft 

enclosure .  40 

39.  Ratio  of  effective  conductivity  to  the  conductivity  of  air  vs.  the 

average  interior  temperature  (4-in.  pipe) .  41 


Go  to  Contents 


iv 


Page 


Figure 

40.  Ratio  of  effective  conductivity  to  the  conductivity  of  air  vs.  the 

average  interior  temperature  (2-in.  pipe) .  41 

41.  Ratio  of  effective  conductivity  to  the  conductivity  of  air  vs.  the 

average  interior  temperature  (4-in.  insulated  pipe) .  42 

42.  Ratio  of  effective  conductivity  to  the  conductivity  of  air  vs.  the 

average  interior  temperature  (2-in.  insulated  pipe) .  42 

43.  Effective  conductivity  correlations .  43 

44.  Correlation  of  pipe  radius  with  intercepts .  43 

45.  Temperature  contours  for  sq4ib65a .  45 

46.  Inside  surface  temperatures  for  sq4ib65a .  45 

47.  Temperature  contours  for  sq2d35c . 45 

48.  Inside  surface  temperatures  for  sq2d35c .  46 

49.  Temperature  contours  for  sq2ic65a .  46 

50.  Inside  surface  temperatures  for  sq2ic65a .  46  * 

51.  Ratio  of  effective  conductivity  to  the  conductivity  of  air  vs.  the 

average  interior  temperature .  47 

52.  Polynomial  curve  fits  to  experimental  boundary  data .  48 

53.  Finite  element  mesh  used  for  conduction  solution  of  the  2-ft  x  4-ft 

enclosure .  48 

54.  Comparison  of  experimental  and  numerical  inside  insulation 

surface  temperatures .  48 

TABLES 

Table 

1.  Methods  of  determining  the  effective  thermal  conductivity  of  an 

air  gap .  12 

2.  Comparison  of  published  velocity  predictions  with  FECOME  for 

a  vertical  enclosure .  22 

3.  Comparison  of  viewfactor  calculations .  24 

4.  Comparison  of  radiation  heat  fluxes .  24 

5.  Test  apparatus  configurations .  26 

6.  Locations  of  thermocouples  in  the  1-ft  x  1-ft  enclosure .  28 

7.  Locations  of  thermocouples  in  the  2-ft  x  4-ft  enclosure .  29 

8.  Enclosure  dimensions  and  effective  gap  widths .  32 

9.  Configurations  for  the  numerical  experiments .  32 

10.  Coefficients  for  eq  133 .  43 

11.  Comparison  of  effective  conductivity  correlations .  44 

12.  Conditions  for  comparisons  between  numerical  conduction  and 

convection  solutions .  44 


v 


Back  to  Contents 


NOMENCLATURE  AND  ABBREVIATIONS 


A 

a 

B 

C 

CP 

Cv 

CRREL 

D 

d.o.f. 

E 

F 

FE 

FECOME 

FERF 

FEVIEW 

G 

S 

Gr 

H 

h 

k 

L 

Le 

N 

n 

Nu 

P 

V 
Pr 

Q 

q 

Qj 

R 

r 

Ra 

^XXX 

s 

T 

t 

U 

V 

w 

X 

V 

y 

z 

a 

P 


A  heat  transfer  correlation  coefficient 
Area 

A  heat  transfer  correlation  coefficient 
A  heat  transfer  correlation  coefficient 
Specific  heat  at  constant  pressure 
Volumetric  specific  heat,  (Cpp) 

Cold  Regions  Research  and  Engineering  Laboratory 
Diameter 

Degrees  of  freedom 
Eccentricity  (eq  50) 

Radiation  viewfactor 
Finite  element 

Finite  Element  Combined  Equations  (computer  code) 

Frost  Effects  Research  Facility 

Computer  code  using  FE  to  determine  radiation  view  factors 
Vertical  gap  width  (eq  48) 

Acceleration  due  to  gravity 

Grashof  number 

Inside  height  of  an  enclosure 

Convective  heat  transfer  conductance 

Thermal  conductivity 

Arc  length  of  an  element  side 

Hypothetical  gap  width  (R-r),  or  characteristic  length 

Interpolation  functions 

number  of  nodes 

Nusselt  number 

Perimeter 

Pressure 

Prandtl  number 

Internal  heat  generation 

Heat  flux 

Radiation  heat  flux  through  surface  j 

Hypothetical  radius  of  a  circle  with  the  same  parameter  as  the  enclosure 

Radius  of  a  cylinder 

Rayleigh  number 

Thermal  resistance  of  xxx 

Boundary  surface 

Temperature 

Thickness 

x-direction  velocity 

y-direction  velocity 

Inside  width  of  an  enclosure 

x-coordinate  location 

Inside  height  of  enclosure 

y-coordinate  location 

z-coordinate  location 

Thermal  diffusivity 

Coefficient  of  thermal  expansion 


Go  to  Contents 


vi 


5kj  Dirac  delta  function,  equal  to  1  if  k=j  and  equal  to  0  if 
Emissivity  of  surface  j 

cp  Heat  flux  across  a  boundary  s,  or  dimensionless  temperature  (eq  30) 
v  Kinematic  viscosity 
\i  Dynamic  viscosity 
p  Density 

Shape  function  parameter  (vertical  direction) 
r|  Shape  function  parameter  (horizontal  direction) 
a  Boltzmann's  constant 


Subscripts 

air 

b 

c 

ci 

cond 

conv 

D 

E 

eff 

eq 

j 

k 

L 

o 

P 

r 

ref 

s 


Of  air 

Boundary  layer  or  distance  traveled  by  the  boundary  layer  on 
cylinder  (nr) 

Convection 

Inner  radius  of  conduit 

Conduction 

Convection 

Diameter 

Exterior  casing 

Effective 

Equivalent 

Indices  for  pipe  number,  directions,  inside  (pipe  diameter  or  radius),  etc. 
Indices  for  pipe  number,  directions,  etc. 

Indices  for  pipe  number,  directions,  etc. 

Characteristic  length,  or  perimeter  lining 
Outside  (diameter  or  radius) 

Pipe 

Radiation 

Reference  for  the  coefficient  of  thermal  expansion 
Sphere 

Reference  for  convective  heat  transfer 


Superscripts 

a  Area 

B  Heat  transfer  correlation  coefficient 
e  Element 

m  Indices  for  iteration  numbers 
p  Pressure 
s  A  boundary  (general) 


vii 


Back  to  Contents 


BLANK 


Go  to  Contents 


Two-Dimensional  Analysis  of 
Natural  Convection  and  Radiation  in  Utilidors 

PAUL  W.  RICHMOND 


INTRODUCTION 

Many  large  building  complexes,  such  as  military  facilities  and  university  cam¬ 
puses,  are  served  by  central  heat  distribution  systems.  Utilidors  are  often  used  to 
contain  the  heat  distribution  lines  and  other  utilities  between  utility  stations  and 
the  serviced  buildings.  These  enclosures  are  generally  constructed  of  concrete  and 
are  usually  installed  below  ground  level.  Other  materials,  such  as  wood  and  sheet 
metal,  are  also  used.  Figure  1  shows  cross  sections  of  two  utilidors  constructed  in 
Arctic  regions. 

In  the  southern  United  States,  utilidors  are  referred  to  as  utility  trenches  and 
often  have  their  upper  side  (lid)  at  ground  level.  Because  utilidors  are  used  to  dis¬ 
tribute  heat  (steam  or  hot  water),  it  is  important  to  know  what  the  heat  loss  from 
the  utilidor  is  in  order  to  estimate  losses  in  the  heat  distribution  system  and  to 
design  for  efficient  use  of  insulation.  Additionally,  the  presence  of  unheated  lines 
(e.g.,  domestic  water,  fire  protection,  or  sewer  lines)  requires  that  the  air  tempera¬ 
ture  within  the  utilidor  remain  above  freezing.  Because  of  the  complexity  of  the 
geometry,  the  heat  transfer  analysis  must  be  done  using  approximate  or  numeri¬ 
cal  procedures. 

Recently,  numerical  methods  (finite  difference,  finite  element)  have  been  used 
to  evaluate  utilidor  performance.  Modeling  of  conductive  heat  transfer  around 
utilidors  and  building  foundations  has  been  done  successfully  (Zarling  and  Braley 
1984,  Phetteplace  et  al.  1986,  Kennedy  et  al.  1988).  Modeling  of  heat  transfer  within 
utilidors  is  generally  done  using  conductive  approximations  (Smith  et  al.  1979). 
Convection  and  radiation  can  have  a  significant  effect  on  the  total  heat  transfer, 
and  accurate  models  are  necessary.  Once  available,  numerical  models  and  corre¬ 
lations  of  heat  transfer  within  utilidors  can  be  used  in  the  design  process  of  new 
utilidors.  A  second  application  is  in  the  thermal  evaluation  of  existing  utilidors 
for  rehabilitation  or  renovation,  replacing  current  approximation  methods. 

Researchers  of  numerical  methods  have  demonstrated  that  convection  and  radi¬ 
ation  can  be  modeled  using  finite  element,  finite  difference,  or  other  numerical 
techniques  (Gebhart  et  al.  1988,  Arpaci  and  Bayazitoglu  1990).  However,  these 
efforts  have  not  been  applied  to  utilidors,  and  in  general  have  been  limited  to  simple 
geometries.  Additionally,  the  two  modes  of  heat  transfer  are  not  often  evaluated 
in  tandem. 
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4-in.  Polyurethane 


as  Required  Heat  Trace  Pipe 

Figure  1.  Cross  sections  of  two  utilidors  constructed  in  the  Arctic.  (After  U.S.  Army  1987.) 


Utilidor  sizes  and  shapes  are  determined  by  considering  the  number  and  sizes 
of  the  pipes  they  will  contain,  their  location  relative  to  the  ground  surface,  and  the 
ease  of  access  desired  for  maintenance  or  repairs.  Phetteplace  et  al.  (1981)  presented 
the  utilidor  and  pipe  sizes  for  all  the  utilidors  located  on  Fort  Wainwright,  Alaska. 
They  reported  approximately  200  different  configurations;  utilidor  sizes  ranged 
from  1  ft  x  1  ft  to  7  ft  x  9  ft,  and  pipe  sizes  varied  from  1  in.  to  24  in.  in  diameter. 
Clearly,  it  is  not  possible  to  conduct  physical  experiments  using  every  combina¬ 
tion  of  utilidor  size  and  pipe  combination. 

The  objective  of  this  work  was  to  investigate  convection  and  radiation  in  enclo¬ 
sures,  specifically  rectangular  utilidors  containing  one  or  more  heated  pipes.  The 
work  presented  considers  the  steady-state,  two-dimensional  problem  of  convec- 
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tion  and  radiation  within  an  enclosure.  Results  of  numerical  and  experimental 
investigations  are  combined  to  obtain  a  methodology  for  the  two-dimensional  ther¬ 
mal  analysis  of  utilidors. 


BACKGROUND 

The  governing  equations  for  incompressible  Newtonian  fluid  flow  in  an  enclo¬ 
sure  are  the  Navier-Stokes  (momentum)  equations,  the  energy  equation,  and  the 
continuity  equation.  The  steady  state,  laminar  flow  momentum  equations  are 


du  du ]  1  dp  f  d2u  d2u 


3t>  dv ^ 

VT~+UT~ 
dy  dxj 


g|5(T-Tref)+ 


1  dp 

pay 


-V 


d2v  d2v ^ 
dx2  +  dy2  > 


=  0 


(2) 


for  a  two-dimensional  flow  field,  where  y  is  the  vertical  direction  and  x  is  the  hori¬ 
zontal  direction.  The  continuity  equation  is 


du  dv 

3^+3y_ 

and  the  energy  equation  (neglecting  viscous  dissipation)  is 


(3) 
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The  energy  equation  reduces  to 


.  (d2T 

[dx2+dy2  J 


-Q  =  o 


(5) 


for  solid  regions  with  homogeneous,  isotropic  materials,  and  constant  thermal  con¬ 
ductivity  ( k ).  These  equations  are  coupled  and  result  in  four  equations  and  four 
unknowns:  pressure,  temperature,  and  the  x  and  y  components  of  velocity  (p,  T,  u, 
and  v).  For  complex  geometries,  these  equations  cannot  be  simplified  and  solved 
directly. 

Heat  transfer  correlations  for  convective  heat  flow  in  enclosures  are  generally 
expressed  in  terms  of  the  Nusselt  number  (Nu)  and  the  Rayleigh  number  (Ra).  These 
dimensionless  parameters  are  defined  as 
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where  g  is  the  acceleration  due  to  gravity,  |5  is  the  thermal  coefficient  of  expan¬ 
sion,  \i  is  dynamic  fluid  viscosity,  v  is  the  kinematic  fluid  viscosity,  hc  is  the  heat 
transfer  conductance,  a  is  the  thermal  diffusivity,  and  k  is  the  thermal  conductiv¬ 
ity  of  the  fluid.  Pr  is  the  Prandtl  number  and  Gr  is  the  Grashof  number.  The  two 
remaining  undefined  terms,  AT  and  L,  are  dependent  on  the  boundary  conditions 
and  geometry  of  the  problem.  In  the  simplest  case,  AT  will  be  the  temperature  dif¬ 
ference  between  a  warm  surface  and  a  cold  surface.  The  variable  L  is  a  character¬ 
istic  length  of  the  geometry.  For  concentric  cylinders  the  difference  in  radii  or  gap 
width  is  often  used;  other  examples  are  discussed  below.  Correlations  for  Nu  are 
found  in  the  form  of 


Nu  =  ARaB 

when  a  specific  material,  such  as  air,  is  specified,  or 


(10) 


Nu  =  AGrB  (H) 

for  the  general  case  of  natural  convection  in  fluids  or  gases. 

Heat  transfer  by  radiation  between  two  surfaces  can  have  a  large  effect  on  the 
heat  transfer  correlations.  Experiments  and  analytical  or  numerical  analysis  can 
include  these  effects  or  they  can  be  removed.  Radiation  is  primarily  reflected  in 
the  heat  transfer  conductance  h.  Generally,  h  should  be  considered  to  be  the  sum 
of  two  components,  hY  and  hc,  i.e.,  the  conductances  due  to  radiation  and  to  con¬ 
vection.  It  is  not  always  clear  when  examining  heat  transfer  correlations  if  this  is 
the  case,  or  if  h  represents  merely  hc. 

A  vertical  rectangular  cavity  (enclosure)  is  defined  as  an  enclosure  bounded  by 
two  vertical  surfaces  held  at  different  temperatures.  The  other  two  parallel  sur¬ 
faces,  top  and  bottom,  are  taken  as  insulated  (Gebhart  et  al.  1988).  Heat  transfer 
occurs  only  at  the  vertical  surfaces.  The  characteristic  length  L  for  this  geometry  is 
the  distance  between  the  hot  and  cold  walls,  and  the  characteristic  temperature 
AT  is  the  difference  between  the  vertical  wall  temperatures.  For  an  air-filled  square 
enclosure,  Ostrach  (1972)  summarized  the  following  numerical  results  for  the 
average  Nusselt  number  in  the  form  of  eq  11. 


Reference 

A 

B 

eq 

Newell  and  Schmidt  (1969) 

0.0547 

0.397 

(12) 

Han  (1967) 

0.0782 

0.3594 

(13) 

Elder  (1965) 

0.231 

0.25 

(14) 

He  reported  tolerable  agreement  between  these  correlations  and  experiments. 

Recently,  de  Vahl  Davis  and  Jones  (1983)  presented  a  numerical  benchmark 
solution  for  air  in  a  square  vertical  enclosure  at  Ra  values  from  103  to  106.  Fitting 
an  equation  to  their  Nusselt  number  data  at  the  cold  surface  yields 


Nu  =  0.141 62Gr0,2996. 


(15) 


Equations  12-15  are  drawn  on  Figure  2. 

Correlations  have  also  been  developed  for  vertical  enclosures  with  aspect  ra¬ 
tios  (height /width)  other  than  one.  Gebhart  et  al.  (1988)  present  several  correla¬ 
tions  of  the  form 
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Figure  2.  Heat  transfer  correlations  for  vertical  enclosures. 


NM  =  AGrBQQC  (16) 

where  Y/W  is  the  height/width  ratio,  and  A,  B,  and  C  are  the  constants  for  air  (see 
table  below). 


Reference 

A 

B 

C 

eq 

Newell  and  Schmidt  (1969) 

0.155 

0.315 

-0.265 

(17) 

Eckert  and  Carlson  (1961) 

0.119 

0.3 

-0.1 

(18) 

Jakob  (1949) 

0.18 

0.25 

-0.111 

(19) 

MacGregor  and  Emery  (1969) 

0.25 

0.25 

-0.25 

(20) 

Horizontal  rectangular  enclosures  are  described  as  cavities  in  which  the  lower 
horizontal  surface  is  heated  while  the  upper  surface  is  cooled;  the  sides  are  insu¬ 
lated.  The  correlations  obtained  by  several  researchers  can  be  presented  in  the  form 
of  eq  10,  when  the  Prandtl  number  for  air  is  taken  as  0.72.  The  characteristic  length 
L  in  the  Ra  number  is  the  height  of  the  enclosure.  The  constants  for  several  corre¬ 
lations  are  shown  in  the  following  table  (Gebhart  et  al.  1988). 

Reference  A  B  eq 

Dropkin  and  Somerscales  (1965)  0.0673  0.3333  (21) 

Silveston  (1958)  0.0877  0.31  (22) 

Kraichnan  (1962)  0.1524  0.3333  (23) 


Probably  the  most  investigated  enclosure  containing  an  interior  heat  source  is 
a  concentric  pipe  system.  Gebhart  et  al.  (1988)  reviewed  the  significant  number  of 
experimental  and  numerical  investigations  of  this  geometry,  noting  that  different 
nondimensional  systems  have  been  used  in  most  studies.  For  correlations  based 
on  mean  heat  transfer  rates,  gap  width  (outer  radius-inner  radius)  is  often  used  as 
the  characteristic  length  (L).  An  example  of  this  is  the  following  equation  by  Grigull 
and  Hauf  (1966): 
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(24) 


Nul  = 


0.2  +  0.145  I  h 


Gr °-25c 


1,02  (%) 


for  30, 000  <GrL<  716, 000 
and  0.55  <  yjj  <  2.65 


where  D{  is  the  diameter  of  the  internal  cylinder.  Gap  width,  however,  does  not 
provide  all  the  heat  transfer  information  that  may  be  desired,  i.e.,  the  conductances 
for  the  two  surfaces  are  not  obtained  individually,  but  are  lumped  together.  The 
results  of  many  studies  are  presented  using  an  equivalent  conductivity  (fceq),  which 
is  defined  as  the  ratio  of  actual  heat  flow  to  that  due  to  conduction  alone  across 
the  region.  For  concentric  cylinders,  the  equivalent  conductivities  based  on  the 
inside  and  outside  surface  areas  are 


eqA  Nucond  2k  [dJ 


Nu0 

Nucond 


Mo  in 
2k 


(25) 

(26) 


where 

Nucond  =  ^(DJDi)  ' 


(27) 


The  total  energy  lost  by  one  cylinder  equals  that  gained  by  the  other  (i.e.,  eq  25 
equals  eq  26).  The  subscript  i  refers  to  the  inner  cylinder  and  o  to  the  outer  one, 
and  Nucond  is  the  Nusselt  number  for  pure  conduction  between  concentric  cylin¬ 
ders  (Gebhart  et  al.  1988). 

Kuehn  and  Goldstein  (1978)  combined  a  large  amount  of  data  and  obtained  the 
following  correlations  for  Pr  =  0.7  (air): 


Nu:=- 


ln  1  +  2  j 


(0.5R««)15+(0.12R»W3) 
-2 


Nu0  =  - 


15' 


1/15 


lm  1-2 


<t»b 


Nu-i  _  (^b  “  ^o) 


1/15 


Nui+Nu0  (Ti-T0) 


^uconv 


VNui  Nuo  J 


Nu, 


cond 


ln(D0/Di) 

Nu  =  [{Nucondf+{Nuconv) 


15 


1/15 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 
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Nu 


(34) 


eCl 


N^cond 


where  the  Nusselt  numbers  are  averaged  values  for  the  overall  heat  transfer  around 
the  cylindrical  surfaces,  and  are  based  on  D0.  RaDi  is  the  Rayleigh  number  based 
on  Di  and  RaDo  is  that  based  on  D0.  The  temperature  difference  in  Ra  is  the  differ¬ 
ence  between  the  inner  (T{)  or  outer  (T0)  surface  temperatures  and  the  average  fluid 
temperature  (Tb)  between  the  inner  and  outer  cylinder  boundary  layers.  Tb  can  be 
determined  from  (f^,  the  average  dimensionless  fluid  temperature  between  bound¬ 
ary  layers.  An  iterative  solution  to  the  correlation  will  be  required  to  obtain  the 
Nusselt  numbers.  What  is  significant  about  this  correlation  is  that  the  conductances 
for  both  surfaces  can  be  obtained  along  with  the  mean  fluid  temperature. 

Lunardini  (1990)  conducted  experiments  using  a  conduit  system  used  at  many 
government  installations  (Fig.  3).  He  identified  four  ways  to  evaluate  the  thermal 
resistance  of  the  air  gap  Ra  given  by 


(35) 


from  the  Federal  Guide  Specification  (1981),  where  the  convective  coefficient  ( h ) 
assumes  a  constant  value  of  3  Btu/hr  ft2oF,  or 

f  \ 


In 


R 

3  2  Jtfceff 


(36) 


where 

fceff  =0.11Ra°'29fcair  (37) 

obtained  from  Grober  et  al.  (1961),  or  from  his  own  data 

fceff  =  1.463 Ra°123fcair/  (38) 

which  includes  radiation  effects,  or 


ke{(  =  0.68Rfl°'157fcair, 

which  has  had  the  effect  of 
radiation  removed.  /ceff  is 
the  effective  conductivity  of 
air,  kair  is  the  conductivity  of 
air,  rci  is  the  inner  radius  of 
the  outer  conduit,  and  rx  is 
the  outer  radius  of  the  insu¬ 
lation.  The  air  gap  thickness 
is  used  as  the  characteristic 
length  in  the  Rayleigh  num¬ 
ber. 

Boyd  (1981)  combined 
data  from  concentric  circu¬ 
lar  cylinders  with  data  from 


(39) 


Carrier  Pipe 
Pipe  Insulation 
Outer  Conduit 


Figure  3.  Cross  section  of  a  concentric  pipe  conduit. 
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hexagonal  cylinders  inside  a  circular  cylinder.  He 
found  that  the  Nu  should  be  based  on  gap  width 
while  Ra  should  be  based  on  the  radius  of  the 
internal  cylinder.  This  approach  indirectly 
includes  the  aspect  ratio  used  by  other  investi¬ 
gators  (e.g.,  eq  24). 

Powe  and  Warrington  (1983)  and  Warrington 
and  Powe  (1985)  investigated  cylinders  and 
spheres  mounted  in  spherical  or  cubical  enclo¬ 
sures.  Although  their  experimental  correlations 
are  probably  not  appropriate  to  this  study  some 
of  their  observations  are  of  interest.  They  used 
parameter  L/rs  as  a  multiplier  to  the  Ra  number 
in  correlations  similar  to  those  above,  where  L  is 
the  gap  width  and  rs  is  the  hypothetical  spheri¬ 
cal  radius  based  on  volume.  This  parameter  is 
used  to  account  for  the  observation  that,  as  the 
interior  body  becomes  smaller,  the  natural  con¬ 
vection  phenomena  can  be  divided  into  three 
regimes.  These  regimes  are  (1)  infinite  atmo¬ 
sphere  solution  for  large  L/rs,  (2)  enclosure  solu¬ 
tions  for  moderate  L/rs,  and  (3)  conduction  solu¬ 
tions  for  small  L/rs.  Additionally,  Warrington  and 
Powe  (1985)  determined  that  for  nonisothermal  figure  4.  Rectangular  enclosure 
internal  bodies,  analyses  using  the  average  body  configuration  ofGhaddar  (1992). 

temperature  compared  well  with  results  from  iso¬ 
thermal  internal  bodies. 

Ghaddar  (1992)  conducted  a  numerical  study  of  a  uniformly  heated  (constant 
heat  flux)  cylinder  in  an  enclosure  as  shown  in  Figure  4.  Note  that  the  pipe  is  not 
centered  vertically,  but  is  in  the  lower  portion  of  the  enclosure.  She  used  a  con¬ 
stant  wall  temperature  of  59°F,  and  varied  the  heat  flux  into  the  cylinder;  a  mean 
cylinder  temperature  was  used  to  calculate  the  Rayleigh  and  Nusselt  numbers.  Her 
numerical  model  did  not  include  radiation.  The  heat  transfer  correlations  devel¬ 
oped  were 

,  n -I  0.207 

Nuh  =  1.81  Rah  —  (40) 

.  VP )_ 

Nub  =  0.604fti°-2083  (41) 

where  L  is  the  hypothetical  gap  width,  rp  is  the  pipe  radius,  and  b  is  the  distance 
traveled  by  the  boundary  layer  on  the  pipe  (half  the  pipe  circumference).  The 
hypothetical  gap  width  is  defined  as  the  difference  between  the  effective  radius  of 
a  cylinder  that  has  a  circumference  equal  to  the  perimeter  of  the  noncircular 
enclosure  and  the  radius  of  the  interior  pipe.  Equation  40  becomes  eq  42  after 
inserting  Ghaddar's  test  conditions  into  the  L/rp  term: 

Nuh=3.756Ral207.  (42) 
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Figure  5.  Experimental  configuration  of  Stewart  and  Verhulst  (1985) 


Stewart  and  Verhulst  (1985)  presented  the  results  of  experiments  in  which  two 
heated  cylinders  were  in  a  cooled  rectangular  enclosure.  Figure  5  shows  their 
apparatus,  which  was  filled  with  distilled  water;  measurements  were  made  with 
both  cylinders  heated  and  when  heated  individually.  They  investigated  a  number 
of  different  characteristic  lengths  and  found  that  the  best  correlation  (least  devia¬ 
tion  from  the  data)  occurred  when  the  hypothetical  gap  width  L  was  used.  (When 
more  than  one  pipe  was  used  to  calculate  L,  an  effective  radius  that  included  both 
interior  pipes  was  used.) 

For  both  cylinders  heated 


Nul  =  0.420Ra°-219 

(L  includes  both  cylinders) 

(43) 

Nuh  =  1.534Rfl°169 

(L  using  large  cylinder  only) 

(44) 

Nul  =  0.231Rfl°' 243 

(L  using  small  cylinder  only). 

(45) 

For  only  one  cylinder  heated 

Nul  =0.256RaJ266 

(large  cylinder  heated,  L  using  large 
cylinder  only) 

(46) 

Nul  =  0.027Ra  £371 

(small  cylinder  heated,  L  using  small 
cylinder  only). 

(47) 

Babus'Haq  et  al.  (1986)  used  interferometric  flow  visualization  to  determine  the 
optimized  location  of  a  single  warm  pipe  in  a  cool  square  enclosure  with  the  antici¬ 
pated  application  being  district  heating  distribution  lines,  i.e.,  utilidors.  Figure  6 
is  a  diagram  of  their  experimental  apparatus.  Although  Babus'Haq  et  al.  did  not 
develop  any  heat  transfer  correlations  per  se,  their  data  for  heat  loss  from  a  cen¬ 
tered  pipe  to  the  enclosure  walls  can  be  represented  by 

Nug  =  0.34Gr°'25  (48) 
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where  the  characteristic  lengths  G  and  D  are  the  average  vertical  gap  width, 
(Y-D)  /2,  and  the  pipe  diameter,  respectively.  This  equation  can  be  converted  to 
the  following  form  using  their  test  conditions: 


NwL=0.4048Ra^25. 


(49) 


Additionally,  they  found  that  the  optimal  location  for  a  heated  pipe  in  a  cooled 
square  enclosure  is  in  the  upper  part  of  the  enclosure,  specifically  at  £  =  -0.73,  where 
E  is  the  eccentricity  given  by 


E  = 


2H  "1  . 

Y-D  J  J 


(50) 


and  H  is  the  distance  from  the  top  of  the  pipe  to  the  inside  of  the  enclosure  lid.  Y  is 
the  interior  vertical  dimension  and  D  is  the  pipe  diameter.  Figure  7  compares  the 


Figure  7.  Heat  transfer  equations  for  pipes  in  enclosures. 
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heat  transfer  correlations  based  on  the  hypothetical  gap  width  L.  All  of  the  equa¬ 
tions  yield  Nusselt  numbers  within  20%  of  each  other  with  the  exception  of 
Ghaddar's  (eq  42),  which  is  about  260%  higher  than  the  mean  value  of  the  other 
equations  at  a  Rayleigh  number  of  108.  This  could  be  due  to  the  pipe  location  ( E  = 
0.52  using  eq  50),  which  agrees  with  the  findings  of  Babus'Haq  et  al.  (1986)  that 
more  heat  transfer  occurs  from  hot  pipes  when  placed  lower  in  the  enclosure  (posi¬ 
tive  values  of  E). 

Currently  accepted  practice  by  Federal  agencies,  for  the  thermal  analysis  of  the 
utilidors  shown  generically  in  Figure  8,  is  presented  by  Smith  et  al.  (1979)  and  by 
the  U.S.  Army  (1987).  Two  assumptions  are  made:  (1)  the  air  temperature  inside 
the  utilidor  is  uniform  and  (2)  interior  air  film  resistance  can  be  ignored.  The  pro¬ 
cedure  consists  of  determining  the  thermal  resistances  by  assuming  that  the  rect¬ 
angular  enclosures  can  be  treated  as  circular  by  using  a  radius  calculated  from  the 
mean  perimeters  (PL  and  PE  in  Fig.  8).  If  the  interior  pipes  are  insulated,  the  con¬ 
duction  resistance  of  the  air  gap  is  neglected.  If  the  interior  pipes  are  uninsulated, 
then  the  resistance  may  be  based  on  both  the  air  film  and  pipe  material.  For  mul¬ 
tiple  pipes  with  differing  temperatures,  all  of  the  resistances  and  pipe  tempera¬ 
tures  are  included  to  obtain  an  interior  air  temperature. 

It  is  also  possible  to  determine  an  effective  conductivity  of  the  air  that  includes 
all  the  film  resistances,  radiation,  and  natural  convection  effects.  These  procedures 
depend  upon  estimates  of  rectangular  enclosures  as  circular  and  neglecting  any 
effects  of  eccentricity  of  the  pipe  location.  These  approaches  are  illustrated  as  fol¬ 
lows:  Using  the  square  enclosure  in  Figure  8,  the  heat  loss  per  unit  length  is 


Q  = 


AT 

1R 


(51) 


where  AT  is  the  difference  between  T0  and  T3  and  SR  is  the  sum  of  the  resistances. 
With  the  assumption  that  the  square  enclosure  can  be  treated  as  a  cylinder  of  equal 
perimeter,  the  resistances  are  determined  as 


*pipe 

“  i,  .  p  .  i 

''pipe  1  pipe  L 


(52) 


Figure  8.  Generic  utilidors  for  current  utilidor  thermal  analysis  procedure. 
(After  U.S.  Army  1987.) 
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P  _  ^insulation 

-''pipe  insulation  t  p  i 

ftincnlaHrm  -*■  in  cilia  firm  '  -t 


(53) 


R*  i, 


—  A  nr  p 
air  gap  -  2jcfceff  •  1 '  °  airgaP 
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A 


where 


2nke{[+A2h[  ln^ 


(54) 


K  = 


°(T12+T2 

)(T1+T2 

pull 

_e2  ^2 

lel  ). 

(55) 


ft 


_  ''lining 
thermal  lining  —  T  p  7 
^lining 1 L  ’ 1 


^exti 


'casing 


enor  casmg  i  pi 

^casing  1 E  1 


(56) 

(57) 


where  k  is  the  conductivity,  P  is  the  mean  perimeter,  D0  is  the  outside  diameter,  Dj 
is  the  inside  diameter  (of  the  air  gap),  and  t  is  the  thickness  of  the  casing  or  lining. 


Table  1.  Methods  of  determining  the  effective  thermal  conductivity  of  an  air  gap. 


Number 


K<# 


Source 


Comments 


0  .llft^k* 
1.463ft! 


eq  37  Based  on  cylinder,  radiation  included* 

eq  38  Based  on  cylinder,  radiation  included1 


D0  In — — 

°-34GCfc--(^y 


4  OAORa^k^ 

5  0.68ft!  "1S7  k^ 

6  keqfcalr 


eq  48 


Holman  (1976) 
eq39 
eq  28-34 


Based  on  rectangular  enclosure,  y  is  enclo¬ 
sure  height,  includes  radiation  ** 

Based  on  cylinder,  radiation  not  included. 

Based  on  cylinder,  radiation  not  included. 

Based  on  cylinder,  radiation  not  included. 


7 


1.81 


Ra, 


p.-Pi 

(  Do 

D„ln  — 

Ui 


eq  40  Based  on  rectangular  enclosure,  radiation 

not  included. 


8 


0.23 


(t  -T  V125 
I0 


V  rF  J 


Smith  et  al.  (1979)  Based  on  cylinder,  pipe  is  uninsulated, 
Ta  is  the  air  temperature,  radiation 
included**.  ke ff  is  zero  if  the  pipes  are  in- 
_ sulated. _ 


*  Emissivities  unknown;  the  correlation  is  based  on  work  reported  in  German,  circa  1930. 

*  Emissivities  were  assumed  to  be  0.5  and  0.9  for  the  insulated  pipe  surface  (two  test  conditions),  0.9 


for  the  enclosure. 

**  Materials  were  copper  pipe  and  polymerized  methyl  methacrylate  for  the  enclosure;  no  surface  treat¬ 
ment  or  level  of  copper  pipe  oxidation  was  reported. 

**  This  is  another  "older"  correlation;  the  underlying  references  were  not  given,  but  may  be  attributed 
to  McAdams  (see  Grober  et  al.  1961,  pp.  320-321). 
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The  effective  thermal  conductivity,  keff,  can  be  determined  using  any  of  the  rela¬ 
tionships  in  Table  1.  In  most  cases  an  iterative  solution  will  be  required  to  deter¬ 
mine  the  air  temperature  upon  which  to  base  thermal  properties,  if  unknown,  and 
the  temperature  of  surface  2.  In  general,  the  air  properties  can  be  evaluated  at  the 
average  interior  surface  temperatures.  If  the  effective  conductivity  relation  includes 
radiation,  then  ht  is  zero  in  eq  54.  For  those  correlations  that  do  not  include  radia¬ 
tion,  appropriate  emissivities  can  be  selected  for  use  in  eq  55;  for  those  that  do, 
information  on  the  emissivities  values  used  to  develop  the  correlations  is  limited; 
the  available  data  are  noted  in  Table  1. 

A  number  of  investigators  (Zirjacks  and  Hwang  1983,  Phetteplace  et  al.  1986, 
Kennedy  et  al.  1988)  measured  temperatures  and  heat  flows  in  and  around 
utilidors.  These  measurements  were  generally  extensions  of  modeling  efforts  and 
analysis  was  limited  to  confirmation  of  the  conduction  models  used  to  predict  soil 
temperatures  around  the  utilidors. 


NUMERICAL  MODEL 

Anumerical  model  using  a  finite-element  approach  was  developed  to  solve  the 
momentum,  energy,  and  continuity  equations  in  two  dimensions  for  the  steady- 
state  case.  The  following  assumptions  were  made: 

•  The  fluid  (air)  is  Newtonian  and  incompressible  within  the  Boussinesq  ap¬ 
proximation.  (Fluid  properties  are  constant,  except  for  density,  which  is  a  func¬ 
tion  of  temperature  and  affects  only  the  buoyancy  term.) 

•  Fluid  flow  is  laminar. 

•  Thermal  conductivity  is  constant  for  each  fluid /material. 

Following  Gartling  (1977)  and  Jaluria  and  Torrance  (1986),  the  momentum  equa¬ 
tions  are,  as  stated  earlier. 


9m  9m 
M9x  +  Z,9y 


1  dp 

+  --L--v 
p  dx 


2m  92m 
-  + 


dx2  9  y2 


=  0 


(58) 


d2v^ 

9y2, 


=  0 


(59) 


where  y  is  in  the  vertical  direction  and  x  is  in  the  horizontal  direction.  The  conti¬ 
nuity  equation  is 

3u  +  3^  =  0  (6°) 

dx  9  y 

and  the  energy  equation  is 


9  T  9T^ 

rd2T  92T"| 

l-M  — 

9  y  dxj 

-k 

<dx2  +  dy2/ 

-Q  =  0, 


(61) 


which  becomes  for  a  solid  region  without  convection 


92t"| 

ay2, 


-Q=o. 


(62) 
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The  dependent  variables  p,  T,  u,  and  v  for  the  general  finite  element  e  are 
approximated  by 


ue  =  'ZNi(x,y)ui  (63) 

i= 1 

Ve  =  'ZNi(x,y)vi  (64) 

f= l 

pe  =  2JNj’(x,y)pi  (65) 

!=1 

Te  =  lNi(x,y)Ti.  (66) 

i-1 


Applying  the  Galerkin  criterion  to  element  e  of  an  m+lth  iterate  of  the  governing 
equations,  the  continuity  equation  becomes  (using  eq  63  and  64  and  dV  =  1  dxdy) 

j  MP.d.Ni. dxdyUi  +  J  N[^^Ldxdyvi  =0  (67) 

^  J  dx  Ae  ay 

where  N}  is  the  transpose  of  By  letting  the  notation  ( a,b )  represent  the  area 
integral  of  ab ,  eq  67  becomes 


m\ 

dx  j 


Ui+W 


pdN, 

»'  3y 


(68) 


A  simplification  is  made  at  this  point  in  that  there  are  no  boundaries  with  pres- 
sure  differences.  Using  integration  by  parts  on  S2  terms  in  eq  58, 59,  and  61  yields 


“■  r2'  r  r’N|W 


aNj  9N,\  l/aNi 


dx'dy  p\dx 


,Nf')pi  -vjN{Vue  •nds  =  0 


m  „  3Ni\  „  /9Ni  3Nj\  /3Nj  3N,\  / 


\  3*  3*  /  \  3y  3y  / 

Ml  Ml\„.  .I/Ml  Mf 


3y/r 


+v\-a7'-a7/Hi +p  \ljf /Pi  # 

-gP(Ni,Nj)T1+gPTref(Ni)  =  0 
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Because  there  will  be  no  forced  convection,  the  velocity  at  the  boundary  sur¬ 
face  s  in  eq  69  and  70  will  be  zero,  thus  these  two  terms  drop  out.  In  the  global 
formulation,  the  equations  representing  velocity  boundary  nodes  will  be  set  to  zero 
and  no  other  velocity  boundary  condition  will  be  allowed. 

The  buoyancy  P  is  defined  by 


^Pref-P" 
-  Tref ) 


(72) 


where  Tref  is  a  reference  temperature  at  which  buoyancy  has  no  effect.  Gebhart 
et  al.  (1988)  suggested  using  the  minimum  boundary  surface  temperature  for  the 
reference  temperature,  and  that  suggestion  was  followed. 

The  (n^ ,Qy  and  the  j^N^n  •  VTe<2s  terms  of  eq  71  represent  heat  generated 
within  an  element  and  the  thermal  boundary  conditions.  For  this  application  it  is 
assumed  that  there  is  no  heat  generated  within  an  element,  thus  this  term  is  elimi¬ 
nated.  Expanding  the  remaining  term  to  account  for  specified  heat  flux  and  con¬ 
vective  boundaries  yields 

J^fcn  •  VTeds  =  ^hN^ds  -  ^hN^N^ds  -  j^N^s  (73) 

where  h  and  are  the  convective  heat  transfer  coefficient  and  associated  tem¬ 
perature  and  <j)  is  the  heat  flux  for  the  boundaries  s. 

Summarizing  the  integrals  required  for  eq  68,  69,  70,  71,  and  73,  the  following 
list  is  obtained: 


N?™i 

>'  dy 

™lnf 

dx  ’  > 


Np 
dy  ’  ) 


/ u™  jq. 

V  ‘  '  y  dx 

Lrn  N.  M 

V I  y  3t / 

l^i  m\ 

\dy  '  dy  / 

IdN^'dNA 

\dx'dx 


IdNj  d_N 
\  dx  '  dy  , 


(74) 

(75) 

(76) 

(77) 

(78) 

(79) 

(80) 

(81) 

(82) 
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IdN,  3Nj\ 

(83) 

\dy  '  dx  / 

(84) 

fa) 

(85) 

JsNi  ds 

(86) 

JsNiNjds. 

(87) 

Interpolation  functions 

In  the  above  equations,  the  matrix  N  represents  interpolation  functions  for  an 
arbitrary  element.  bF  are  the  interpolation  functions  one  order  lower  than  N.  Any 
two-dimensional  shape  element  can  be  used  for  this  set  of  equations  so  long  as  C° 
continuity  is  maintained 
(Huebner  and  Thornton 
1982).  Considering  that  most 
of  the  utilidor  components 
are  rectangular  in  shape 
(walls  and  insulation),  it 
would  be  convenient  to  use 
rectangular  elements.  How¬ 
ever,  the  presence  of  pipes 
requires,  at  a  minimum,  tri¬ 
angular  elements  to  model 
these  curved  surfaces.  Many 
triangular  elements  will  be 
required  to  model  the  pipes 
and  meld  the  curved  areas  to 
rectangular  areas.  By  using  a 
rectangular  element,  which  can  have  curved  sides,  fewer  elements  will  be  required. 

Ergatoudis  et  al.  (1968)  presented  the  interpolation  functions  for  curved 
isoparametric,  quadrilateral  elements.  An  element  of  this  shape  is  shown  in  Fig¬ 
ure  9;  the  element  is  defined  by  eight  nodes,  three  on  each  side.  The  interpolation 
functions  are 

for  nodes  1, 3,  5,  and  7 

S=±i,  ii=±i  <88) 

for  nodes  4  and  8 

Ni(tii)=^(i+^2)(i+'n'ni)  £=o,  t)=±i  (89) 

and  for  nodes  2  and  6 

Ni(^Ti)=|(1+^i)(1_'n2)  ^=±i,ti=o.  (9°) 


5 


Figure  9.  Curved  isoparametric  quadrilateral  element. 
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The  variables  2;  and  T)  are  local  variables;  for  an  individual  element  they  are 
related  to  the  global  xi  and  y{  coordinates  by 


x=YlNi  (91) 

i= 1 

y=SNi(^ii)yi.  (92) 

i= 1 


The  derivatives  3N/32,  and  3N/3r|  can  be  found  directly;  however,  these  deriva¬ 
tives  must  be  related  to  £  and  r|  in  order  to  integrate  eq  74  through  87.  This  is  done 
using  the  chain  rule  of  differentiation  and  the  Jacobian  matrix  J;  the  following  rela¬ 
tionship  is  obtained  (Huebner  and  Thornton  1982): 

'dNj  dNj 
dx  3^ 

dNj  dNj 
dy  3il 

Using  the  relationship  dxdy  =  det  Jd  the  above  area  integrals  can  all  be  writ¬ 
ten  in  terms  of  £  and  r|  and  integrated  from  -1  to  1  using  Gaussian  quadrature. 

The  interpolation  functions  (Np)  must  be  linear  (one  order  lower  than  N).  The 
same  element  as  in  Figure  9  is  used;  however,  the  sides  are  assumed  to  be  straight 
and  the  element  is  defined  only  by  nodes  1,  3,  5,  and  7.  The  interpolation  func¬ 
tions  are 

Ni  =-|(1+^i)(1+imi)  (94) 

where  £,  and  t|  take  oh  their  nodal  values  (Fig.  9  and  eq  88-90).  The  evaluation  of 
the  derivatives  and  integrals  follows  the  same  procedure  as  above. 

The  surface  integral,  eq  86,  must  also  be  expressed  in  terms  of  the  parametric 
variables  ^  and  %  and  the  integration  carried  out  over  the  boundary  specified.  In 
order  to  simplify  programming  it  is  assumed  here  that  the  boundary  s  is  made  up 
of  at  least  one  full  side  of  an  element;  thus  from  Figure  9,  side  1  is  described  by 
nodes  1, 2,  and  3;  side  two  is  nodes  3, 4, 5;  side  three  is  nodes  5, 6, 7;  and  side  4  is 
nodes  7,  8,  and  1.  In  this  development  no  other  combinations  are  allowed;  how¬ 
ever,  more  than  one  side  per  element  can  be  specified  as  a  boundary  segment.  For 
each  side  either  %  and  t|  will  be  a  constant  and  ds  is 

ds  =  ~Ledf\  or  ds  =  ~Led£,  (95) 

where  Le  is  the  length  of  the  side.  The  integral  is  now  evaluated  from  -1  to  1,  using 
Gaussian  quadrature.  The  integration  of  eq  87  is  carried  out  similarly,  except  that 
the  term  Nj  is  a  two-dimensional  matrix. 

Solution  procedure 

The  computer  model  FECOME  (Finite  Element  COMbined  Equations,  Rich¬ 
mond  1995)  solves  eq  68-71  simultaneously  for  u,  v,  T,  and  p  and  uses  either  direct 
substitution  or  the  Newton-Raphson  iteration  procedure.  The  solution  procedure 
requires  the  use  of  a  previous  solution  (the  old  solution)  or  an  initial  estimate,  which 
is  then  used  to  obtain  a  new  solution.  Between  iterations,  both  the  direct  substitu- 


=  i  =  U  2, ...  8. 


(93) 
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tion  and  the  Newton-Raphson  method  can  utilize  a  relaxation  procedure,  which 
consists  of  determining  the  weighted  average  of  the  old  and  new  solutions.  The 
equation  is 

new  solution  =  Q(old  solution)  +  (1  -  0)  new  solution  (96) 

where  the  weighting  or  relaxation  factor  (0)  varies  between  0.005  and  0.25  depend¬ 
ing  on  the  maximum  amount  of  change  from  the  previous  solution  and  the  solu¬ 
tion  method.  This  range  was  determined  by  trial  and  error  in  an  effort  to  improve 
the  convergence  rate.  No  formal  optimization  approach  was  attempted,  and  these 
values  are  not  necessarily  the  best  values.  High  values  caused  oscillations  in  the 
direct  substitution  solutions  to  high  Rayleigh  number  problems  (greater  than  105) 
for  the  vertical  enclosure  problem,  and  once  large  oscillations  begin,  the  procedure 
will  not  converge  to  a  solution. 

The  procedure  was  considered  to  have  converged  to  the  steady-state  solution 
when  the  largest  change  in  each  variable  between  successive  solutions  was  less 
than  0.01%.  Changes  this  small  or  smaller  were  found  to  produce  no  significant 
difference  in  the  heat  flux  calculations  through  the  enclosure  sides. 

The  global  matrix  is  3 n+p  by  3n+p  for  the  fluid  elements  plus  n  by  n  for  the  ele¬ 
ments  that  are  a  solid  material,  where  n  is  the  number  of  nodes  and  p  is  the  num¬ 
ber  of  nodes  associated  with  the  pressure  formulation  (four  per  element) .  This  cal¬ 
culation  of  matrix  size  is  reduced  by  the  number  of  solid-fluid  boundary  nodes 
(which  were  counted  twice  in  the  above  analysis).  There  are  28  degrees  of  free¬ 
dom  for  each  element  specified  as  a  fluid  and  8  degrees  of  freedom  for  those  speci¬ 
fied  as  a  solid. 

The  global  matrix  for  the  direct  substitution  method  has  the  form 


~AA  0  0  0 

0  A3  A8  A4 
A1  A9  A7  A5 
[  0  A4t  A5t  0 


'T 

R2 

u 

0 

V 

R1 

JP. 

0 

where 

AA  =  k 

+  C, 


/dNj  9Nj\  /9N|_  <^V 
\ 3*  '  3*  /  \ dy  ’  dy , 


(97) 


(98) 


(99) 


(100) 
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(101) 


A4  =  - 


HlT'N! 

p  \  dx  > 


A5  =  - 


1 

P 


A4T  =- 


1 

P 


(102) 

(103) 


(104) 


A8  =  i) 


dNj  9Nj\ 
dy  '  dx  j 


A9  =  u 


dNj  3Nj\ 
dx  '  dy  / 


Rl  =  ^Tref(Ni) 
Al-tf^Nj) 


R2  =  js  hN^ds  -  Js  <()Nj  ds . 

The  Newton-Raphson  method,  in  its  general  one-dimensional  form,  is 


(D  =  CO0 


/(cop) 

/'(«>o) 


(105) 

(106) 

(107) 

(108) 

(109) 

(110) 


where  (O  is  the  root  of  the  function /(Hombeck  1975).  In  multidimensional  form, 
following  Gartling  (1987), 

new  solution  =  old  solution — /-1  (old  solution)R(old  solution)  (111) 

where  /_1  is  the  inverse  of  the  Jacobian  matrix  of  eq  68-71  and  R  is  the  vector  of  the 
residuals  obtained  by  substituting  the  old  solution  into  eq  68—71.  The  Jacobian 
matrix  is 


3  Rx 

3  Rrp 

3Rt 

o 

dT 

d  u 

dv 

0 

dRu 

dRu 

dRu 

du 

dv 

dp 

dRv 

d  Rv 

dRv 

dRv 

dT 

du 

dv 

dp 

0 

dRp 

dRp 

0 

3  u 

dv 

where  RT ,  Ru,  Rv,  and  Rp  are  eq  68-71,  respectively. 


(112) 
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The  material  properties  of  the  fluid  (air)  can  be  held  constant  or  reevaluated 
between  iterations  using  an  average  temperature  obtained  using  a  number  of 
schemes.  FECOME  averages  the  temperatures  of  the  zero  velocity  nodes  (the 
inside  surfaces  of  the  enclosure)  and  calculates  new  air  properties  based  on  this 
average  temperature  between  each  iteration. 


Radiation 

Large  temperature  differences  can  sometimes  exist  between  utilidor  steam  lines 
and  the  utilidor  walls.  Temperature  differences  between  surfaces  cause  heat  flow 
via  radiation  in  addition  to  natural  convection.  The  heat  flow  due  to  radiation 
(radiosity)  between  surfaces  is  described  by  this  equation: 


h  i 


/§kj  R 

~7. 

VeJ 


n 


y=i 


(113) 


where  a=  Boltzmann's  constant, 

T  =  absolute  temperature  of  surface  k  or  j, 

FM  =  viewfactor  of  surface  k  to  j, 

Qj  =  radiation  heat  flux  into  or  out  of  surface  ;, 
Uj  =  area  of  surface  j , 

£j  =  emissivity  of  surface  j, 

8kj  =  1  if  k  =  ;,  and  =  0  if  k 


The  calculation  of  the  radiation  heat  flux  requires  the  cal¬ 
culation  of  the  radiation  viewfactors  between  each  radia¬ 
tion  surface.  There  are  a  number  of  procedures  to  make 
these  calculations  (Siegel  and  Howell  1992).  Emery  et  al. 

(1991)  made  accuracy  comparisons  between  several  numeri¬ 
cal  approaches.  However,  none  of  the  procedures  are  trivial 
for  complex  geometries.  For  the  two-dimensional  analysis 
of  utilidors,  the  surfaces  should  be  considered  infinite  in 
depth.  By  using  the  finite  element  boundaries  as  the  edges 
of  infinite  strips,  a  special  case  of  two-dimensional  geom-  -  6  — * 

etry  is  obtained.  Figure  10.  Viewfactor 

A  relatively  simple  method  can  be  used  to  obtam  the  analysis  of  F 
viewfactors  for  the  case  of  infinite  strips;  known  as  Hottel's 
crossed-string  method  (Siegel  and  Howell  1992),  the  procedure  is  developed  as 
follows.  To  obtain  the  viewfactor  between  surfaces  1  and  6  in  Figure  10,  first  form 
the  triangle  abc  with  the  infinite  strips  1, 2,  and  3.  The  viewfactors  between  these 
three  surfaces  can  be  written  as 


h-2  +  h-3  =  1 
F2_i  +  F2_3  =  1 
F3-1+F3-2  =  1- 

Multiply  each  equation  by  the  area  of  its  surface: 


aih-2 + aih-3  =  fli 


«2Fj-l  +fl2^>-3  =a2 


(114) 

(115) 

(116) 

(117) 

(118) 
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fl3^3-l  +  fl3^3-2  ~  fl3- 


(119) 


Substituting  the  reciprocity  relations, 

a2^2-l  =  ^1^1—2  (120) 

a3^3-l  =  <Mtl-3  0-21) 


and  solving  the  three  equations  for  F1_2  yields 

_  ai  +  a2~  a3  ,  (122) 

1-2  2«! 

Similarly  for  the  triangle  adb 

p._.  =  al  +g4  ~ a5  .  (123) 

2^q 

Noting  that 

^1-2  +  ^1-4  +  ^1-6  =  1  (124) 

and  solving  eq  122-124  for  F 1-6  yields 

_  +  fl5  ~  fl3  ~  a4  .  (125) 

1-6  2ax 

This  procedure  is  implemented  in  the  program  FEVIEW  (Richmond  1995);  also 
included  is  a  routine  to  check  for  the  shadowing  of  surfaces.  A  surface  is  consid¬ 
ered  shadowed  if  a  line  connecting  the  midpoints  of  two  surfaces  is  intersected  by 
another  radiation  surface.  No  effort  is  made  to  distinguish  partially  shadowed 
elements,  and  as  long  as  the  midpoints  can  be  connected  without  interference,  the 
viewf actor  is  calculated  using  HotteTs  method.  The  viewf actors  are  obtained  prior 
to  running  FECOME  and  appended  to  the  FECOME  grid  data  file.  AFECOME  sub¬ 
routine  uses  eq  113,  nodal  temperatures  and  the  viewfactors,  to  obtain  the  radia¬ 
tion  heat  flux  into  or  out  of  each  of  the  radiation  surfaces. 

The  radiation  heat  fluxes  are  recalculated  at  each  iteration  in  FECOME  using 
the  average  nodal  temperatures  for  each  surface  specified  as  a  radiation  bound¬ 
ary.  In  the  global  formulation,  the  radiation  flux  is  handled  in  the  same  manner  as 
a  boundary  heat  flux  (<|>)  in  eq  73. 


Model  verification 

Verification  of  the  model  consisted  of  comparing  the  model  output  to  known 
(analytical)  or  benchmark  numerical  solutions.  Three  types  of  verifications  were 
done  to  confirm  that  the  model  was  producing  accurate  results;  these  are  described 
in  the  following  paragraphs. 

Several  computer  runs  were  made  to  verify  the  energy  equation  alone  and  the 
implementation  of  the  thermal  boundary  conditions.  These  runs  also  served  to  test 
the  matrix  assembly  and  inversion  routines.  First,  a  square  grid  was  constructed  in 
which  all  the  elements  were  specified  as  a  solid  material  and  two  opposite  sides 
were  set  at  different  temperatures,  with  the  other  two  sides  having  unspecified 
boundary  conditions  (this  corresponds  to  a  zero  heat  flux  boundary).  An  exact 
solution  to  this  simple  one-dimensional  problem  was  obtained.  A  second  test  in  this 
phase  was  a  two-dimensional  conduction  problem;  here  two  adjacent  sides  were 
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set  at  a  constant  temperature 
boundary  and  as  a  thermal 
convection  boundary,  with 
the  remaining  sides  having  a 
zero  heat  flux.  The  results  of 
this  test,  compared  with  the 
analytical  solution  given  by 
Ozisik  (1980),  are  shown  in 
Figure  11;  good  agreement 
was  achieved.  A  third  set  of 
tests  was  run  to  confirm  the 
correct  implementation  of  the 
heat  flux  boundary  condition. 

This  was  done  by  modeling  a 
square  solid  material  with  one 
side  at  a  constant  tempera¬ 
ture,  the  opposite  side  with  a 
specified  heat  flux,  and  the  Fi$ure  lh  Two-dimensional  conduction  problem. 
remaining  sides  unspecified  Dashed  lines  rePresent  ^  analytical  solution,  solid  lines 

/  i.  Tn  x  -m  .  r-  represent  the  numerical  solution. 

(zero  heat  flux).  This  configu-  r 

ration  was  repeated  so  that  all  four  directions  were  tested  with  both  positive  (out 
of  an  element)  and  negative  (into  an  element)  heat  flows. 

In  order  to  test  the  solution  of  the  momentum  equations  and  the  continuity  equa¬ 
tion,  and  their  interaction  with  the  energy  equation,  a  comparison  with  a  well- 
documented  benchmark  numerical  solution  was  done.  As  mentioned  earlier, 
de  Vahl  Davis  and  Jones  (1983)  presented  benchmark  solutions  for  natural  con¬ 
vection  in  a  vertical  enclosure  and  compared  their  results  with  other  investiga¬ 
tors.  A  vertical  enclosure  is  a  closed  cavity  in  which  the  horizontal  surfaces  are 
insulated  (zero  heat  flux  boundaries)  and  the  vertical  sides  are  held  at  two  differ¬ 
ent  temperatures.  His  solutions  were  for  air  at  Rayleigh  numbers  of  103,  104,  105, 
and  106.  Table  2  compares  the  velocity  maximums  along  the  x  =  0.5  and  y  =  0.5 
locations.  Good  agreement  is  observed,  with  percent  differences  less  than  1%.  Also 


Distance  (ft) 
Boundary  at  20°  F 


Table  2.  Comparison  of  published  velocity  predictions  with 
FECOME  for  a  vertical  enclosure. 


Source 

Maximum  velocity 
@x  =  0.5 

Maximum  velocity 
@  y  =  0.5 

y  coordinate 

x  velocity 

x  coordinate 

y  velocity 

Benchmark*  Ra=l(P 

0.813 

3.649 

0.178 

3.697 

Gartling* 

0,824 

3.696 

0.176 

3.696 

FECOME 

0.825 

3.640 

0.175 

3.697 

Benchmark  Ra  =  104 

0.823 

16.178 

0.119 

19.617 

Gartling 

0.824 

16.186 

0.119 

19.630 

FECOME 

0.825 

16.185 

0.125 

19.601 

Benchmark  Ra  - 105 

0.855 

34.73 

0.066 

68.59 

Gartling 

0.854 

34.74 

0.068 

68.63 

FECOME 

0.850 

34.71 

0.067 

68.63 

Benchmark  Ra  =  106 

0.850 

64.63 

0.038 

219.36 

Gartling 

0.854 

64.37 

0.043 

218.43 

FECOME 

0.850 

64.76 

0.033 

218.58 

*  From  de  Vahl  Davis  and  Jones  1983. 
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a .  Temperature  contours.  b.  Stream  lines. 

Figure  12.  FECOME  results  for  a  vertical  square  enclosure  at  Rayleigh  number  105. 


Figure  13.  Grid  used  for  comparing  viewfactors  and  radiation  heat  flux  calcu¬ 
lations. 


shown  are  the  values  submitted  by  Gartling  (to  de  Vahl  Davis  and  Jones).  Gartling 
used  a  model  similar  to  FECOME,  but  with  a  different  mesh.  Figure  12  shows  the 
isotherms  and  streamlines  for  a  Rayleigh  number  of  105;  these  agree  well  with  those 
presented  by  de  Vahl  Davis  and  Jones  (1983). 

Radiation  viewfactor  calculations  obtained  using  the  computer  program 
FEVIEW  were  checked  by  constructing  a  relatively  simple  mesh  (Fig.  13)  and 
using  viewfactor  algebra  and  tables  of  viewfactor  equations  (Siegel  and  Howell 
1992)  to  determine  the  viewfactors  both  manually  and  using  FEVIEW.  Results  of 
this  analysis  are  shown  in  Table  3.  Good  agreement  was  obtained. 

Implementation  of  the  radiation  heat  flux  boundaries  was  checked  by  solving 
eq  113  for  the  heat  fluxes  through  each  of  the  surfaces  in  Figure  13.  A  shell  around 
the  FECOME  subroutine  RADIATE  was  used  for  the  input/output  requirements. 
The  manual  solution  (the  solution  of  the  simultaneous  equations  was  obtained 
using  a  spreadsheet)  is  compared  with  the  computer  solution  in  Table  4;  good  agree¬ 
ment  can  be  seen. 
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Table  3.  Comparison  of  view- 
factor  calculations. 

Surface  1 


to  surface: 

FEVIEW 

Algebraic 

1 

0.000000 

0.000000 

2 

0.000000 

0.000000 

3 

0.000000 

0.000000 

4 

0.000000 

0.000000 

5 

0.000000 

0.000000 

6 

0.000000 

0.000000 

7 

0.000000 

0.000000 

8 

0.000000 

0.000000 

9 

0.004405 

0.004404 

10 

0.012544 

0.012548 

11 

0.018935 

0.018936 

12 

0.023108 

0.023104 

13 

0.015429 

0.015427 

14 

0.021588 

0.021589 

15 

0.030854 

0.030855 

16 

0.044708 

0.044707 

17 

0.064495 

0.064496 

18 

0.089417 

0.089417 

19 

0.112962 

0.112962 

20 

0.123106 

0.123106 

21 

0.019586 

0.019586 

22 

0.036895 

0.036893 

23 

0.089073 

0.089075 

24 

0.292893 

0.292893 

Sum 

0.999998 

0.999998 

EXPERIMENTAL  PROCEDURE 

Experimental  apparatuses 

Two  10-ft-long  experimental  appara¬ 
tuses  were  constructed  to  simulate  sec¬ 
tions  of  typical  rectangular  utilidors.  The 
first  had  an  internal  1-ft  x  1-ft  square  cross 
section;  the  second  was  2  ft  x  4  ft.  Figure 
14  is  the  square  cross  section  apparatus 
with  its  lid  off,  with  a  4-in.  nominal  diam¬ 
eter  pipe  installed.  Heat  transfer  panels 
surrounded  the  sides  of  the  enclosures 
and  a  coolant  was  pumped  through  them. 
The  interior  pipe  was  filled  with  high 
temperature  hydraulic  oil  and  heated  by 
an  internal,  10-ft-long,  1-kilowatt  heating 
element.  The  interior  of  the  enclosure  was 
lined  with  plywood  and  expanded  poly¬ 
styrene  (EPS)  insulation;  the  conductivity 
of  two  samples  of  the  EPS  was  measured 
according  to  ASTM  standards  using  a 
Rapid  K  apparatus.  The  results  of  these 
tests  are  plotted  in  Figure  15.  Thermo¬ 
couples  were  placed  on  either  side  of  the 
EPS  when  installed  in  the  apparatus.  The 


Table  4.  Comparison  of  radi¬ 
ation  heat  fluxes. 


Heat  flux,  Btu/hr 


Surface  no. 

Algebraic 

Numerical 

1 

49.33 

49.33 

2 

64.58 

64.58 

3 

85.71 

85.71 

4 

112.47 

112.47 

5 

-112.47 

-112.47 

6 

-85.71 

-85.71 

7 

-64.58 

-64.58 

8 

-49.33 

-49.33 

9 

-98.66 

-98.66 

10 

-111.42 

-111.42 

11 

-111.42 

-111.42 

12 

-99.32 

-99.32 

13 

-49.33 

-49.33 

14 

-64.58 

-64.58 

15 

-85.71 

-85.71 

16 

-112.47 

-112.47 

17 

112.47 

112.47 

18 

85.71 

85.71 

19 

64.58 

64.58 

20 

49.33 

49.33 

21 

99.32 

99.32 

22 

111.42 

111.42 

23 

111.42 

111.42 

24 

99.32 

99.32 

Sum 

0.66 

0.66 

Figure  14. 1-ft  xl-ft  experimental  appa¬ 
ratus  without  lid. 
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Figure  15.  Thermal  conductivity  of  expanded  polystyrene. 


Plywood  (0.5  in.)  (all  4  sides) 

Figure  16.  Schematic  cross  section  of  the  1-ft  xl-ft  apparatus. 


exterior  of  the  apparatus  was  insulated  and  covered  with  plywood.  Figure  16  is  a 
cross  section  of  the  design. 

Initial  tests  were  done  in  CRREL's  Frost  Effects  Research  Facility  (FERF),  which 
supplied  a  glycol  solution  as  cold  as  -22°F  to  the  apparatus.  During  these  initial 
tests,  it  was  determined  that  significantly  colder  coolant  temperatures  would  be 
required  to  obtain  temperatures  typical  of  Alaska  design  conditions  (-65°F  for  out¬ 
door  air  temperature)  and  to  obtain  near-freezing  temperatures  within  the  enclo- 
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Mean  Temperature  (°F) 

Figure  17.  Thermal  conductivity  of  fiberglass  pipe  insulation. 


Table  5.  Test  apparatus  configurations. 


Enclosure  size 

Nominal  pipe  diam . 

Pipe  treatment 

1  ft  X  1  ft 

4  in. 

uninsulated,  painted* ,  insulated 

2 

uninsulated,  insulated 

2  ft  x  4  ft 

8,4 

insulated 

*  The  uninsulated  pipe  was  painted  with  a  low-emissivity  paint  (Rust- 
Oleum  Aluminum  No.  7715). 


sure.  Unfortunately,  CRREL  had  at  this  time  recently  done  away  with  its  extreme 
low  temperature  capability  for  environmental  reasons. 

Several  years  after  these  initial  experiments,  CRREL  regained  its  extreme  low 
temperature  brine  capability  and  the  apparatus  was  moved  and  replumbed  to  take 
advantage  of  a  coolant  as  low  as  -70°F.  Prior  to  resuming  experiments,  the  inte¬ 
rior  surface-mounted  thermocouples  were  replaced  with  30-gage  surface-mount 
thermocouples  in  order  to  measure  the  surface  temperatures  more  accurately 

Once  experiments  were  resumed,  further  experiments  were  conducted  using 
the  uninsulated  4-in.  pipe.  When  these  were  completed,  experiments  continued 
with  various  pipe  treatments  and  configurations.  Figure  17  is  a  plot  of  the  thermal 
conductivity  of  the  fiberglass  pipe  insulation  used.  Table  5  summarizes  the  test 
apparatus  configurations. 

The  second  apparatus  was  constructed  similarly;  however,  no  plywood  sepa¬ 
rated  the  cooling  panels  from  the  insulation,  and  a  metal  interior  frame  was  used 
to  help  support  the  cooling  panels.  Figure  18  shows  the  2-ft  x  4-ft  apparatus  prior 
to  installing  the  lid.  Figure  19  shows  a  cross  section  with  dimensions  and  pipe 
locations. 

Data  acquisition  system 

Two  different  data  acquisition  systems  were  used,  one  for  the  1-ft  x  1-ft  enclo¬ 
sure  and  another  for  the  2-ft  x  4-ft  enclosure;  type-T  thermocouples  were  used  to 
measure  temperature,  and  a  power  meter  was  used  to  measure  the  power  sup¬ 
plied  to  the  pipes. 

For  the  1-ft  x  1-ft  enclosure,  a  personal-computer-  (PC-)  based  data  acquisition 
system  was  assembled  using  an  80286  processor-based  computer  in  conjunction 
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Figure  18.  2-ft  x4-ft  enclosure  with  4-in.  and  8-in.  insulated  pipes. 


with  Industrial  Computer  Source  data  acquisition  boards.  These  boards  (a  total  of 
four)  were  mounted  in  a  separate  enclosure  and  were  accessed  using  an  8-channel 
multiplexor  board  mounted  in  one  of  the  PC  slots.  Figure  20  is  a  schematic  dia¬ 
gram  of  the  data  acquisition  system;  note  that  an  electronic  ice  point  bath  was  in¬ 
cluded  in  the  thermocouple  circuit  for  temperature  compensation.  The  ice  point 
bath  was  added  because,  during  calibration  of  the  system,  it  was  found  that  the 
onboard  electronic  temperature  compensators  were  not  accurate.  A  data  acquisi¬ 
tion  and  display  program  was  written  to  display  and  store  the  thermocouple  data 
and  measurements  of  the  energy  input  to  the  pipe  heater(s).  These  data  were  stored 
in  ASCII  format  on  floppy  disks  for  later  analysis.  Thermocouple  locations  for  the 
1-ft  x  1-ft  enclosure  are  in  Table  6. 

The  data  acquisition  system  for  the  2-ft  x  4-ft  apparatus  was  based  on  a  Campbell 
Scientific  CR-10  system  with  four  multiplexor  expansion  boards.  One  hundred  and 
twenty-two  thermocouples  were  installed.  Table  7  contains  the  thermocouples'  x 
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Table  6.  Locations  of  thermocouples  in  the  1-ft  x  1-ft  enclosure. 


Channel  no . 

x  coord* 
(in.) 

y  coord* 
(in.) 

Description  of  location* 

0 

3  o'clock  on  pipe 

1 

12  o'clock  on  pipe 

2 

6  o'clock  on  pipe 

3 

9  o'clock  on  pipe 

4 

1.5 

-3.5 

in  air  at  4.5  ft 

5 

6.5 

-3.5 

in  air  at  4.5  ft 

6 

11.5 

-3.5 

in  air  at  4.5  ft 

7 

1.5 

3.5 

in  air  at  4.5  ft 

8 

6.5 

3.5 

in  air  at  4.5  ft 

9 

11.5 

3.5 

in  air  at  4.5  ft 

10 

cold  coolant  inlet 

11 

right  return 

12 

left  return 

13 

lid  return 

14 

bottom  return 

15 

12  o'clock  on  pipe  at  2  ft 

16 

-6.0 

11.0 

outside  insulation  on  left  side 

17 

-5.5 

11.0 

inside  insulation  on  left  side 

18 

6.0 

2.0 

outside  insulation  on  right  side 

19 

5.5 

2.0 

inside  insulation  on  right  side 

20 

6.0 

5.0 

outside  insulation  on  right  side 

21 

5.5 

5.0 

inside  insulation  on  right  side 

22 

6.0 

8.0 

outside  insulation  on  right  side 

23 

5.5 

8.0 

inside  insulation  on  right  side 

24 

6.0 

11.0 

outside  insulation  on  right  side 

25 

5.5 

11.0 

inside  insulation  on  right  side 

26 

6.0 

6.5 

outside  insulation  on  right  side  at  6  ft 

27** 

5.5 

6.5 

inside  insulation  on  right  side  at  6  ft 

28 

-6.0 

6.5 

outside  insulation  on  left  side  at  6  ft 

29** 

-5.5 

6.5 

inside  insulation  on  left  side  at  6  ft 

30 

0.0 

0.0 

outside  insulation  on  bottom  at  6  ft 

31** 

0.0 

0.5 

inside  insulation  on  bottom  at  6  ft 

32 

0.0 

0.0 

outside  insulation  on  bottom 

33 

0.0 

0.5 

inside  insulation  on  bottom 

34 

3.0 

0.0 

outside  insulation  on  bottom 

35 

3.0 

0.5 

inside  insulation  on  bottom 

36 

5.5 

0.0 

outside  insulation  on  bottom 

37 

5.5 

0.5 

inside  insulation  on  bottom 

38 

-3.0 

0.0 

outside  insulation  on  bottom 

39 

-3.0 

0.5 

inside  insulation  on  bottom 

40 

-5.5 

0.0 

outside  insulation  on  bottom 

41 

-5.5 

0.5 

inside  insulation  on  bottom 

42 

-6.0 

2.0 

outside  insulation  on  left  side 

43 

-5.5 

2.0 

inside  insulation  on  left  side 

44 

-6.0 

5.0 

outside  insulation  on  left  side 

45 

-5.5 

5.0 

inside  insulation  on  left  side 

46 

-6.0 

8.0 

outside  insulation  on  left  side 

47 

-5.5 

8.0 

inside  insulation  on  left  side 

48 

0.0 

12.0 

outside  insulation  on  lid 

49 

0.0 

11.5 

inside  insulation  on  lid 

50 

3.0 

12.0 

outside  insulation  on  lid 

51 

3.0 

11.5 

inside  insulation  on  lid 

52 

5.5 

12.0 

outside  insulation  on  lid 

53 

5.5 

11.5 

inside  insulation  on  lid 

54 

-3.0 

12.0 

outside  insulation  on  lid 

55 

-3.0 

11.5 

inside  insulation  on  lid 

56 

-5.5 

12.0 

outside  insulation  on  lid 

57 

-5.5 

11.5 

inside  insulation  on  lid 

58 

0.0 

12.0 

outside  insulation  on  lid  at  6  ft 

59** 

0.0 

11.5 

inside  insulation  on  lid  at  6  ft 

60 

inside  insulation  at  tail 

61 

outside  insulation  at  tail 

62 

inside  insulation  at  head 

63 

outside  insulation  at  head 

*  The  x  origin  is  at  the  center  of  the  enclosure;  the  y  origin  is  at  the  outside  of 
the  bottom  insulation  panel. 

+  Distances  are  measured  from  the  end  that  had  inlets  and  outlets  for  the  cool¬ 
ing  panels;  this  end  is  designated  the  "head"  and  the  opposite  end,  the  "tail." 
Unless  otherwise  stated,  the  thermocouples  were  located  5  ft  from  the  head. 

**  These  thermocouples  were  moved  to  3, 9, 6,  and  12  o'clock  positions,  respec¬ 
tively,  on  the  pipe  insulation  when  it  was  installed;  they  were  not  attached  to 
the  2-in.  bare  pipe. 
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Table  7.  Locations  of  thermocouples  in  the  2-ft  x  4-ft  enclosure. 


T.C.  no. 

x  coord,  (in.) 

y  coord,  (in.) 

T.C.  no. 

x  coord,  (in.) 

y  coord,  (in.) 

Inside  of  bottom  insulation  panel 

61 

24.5 

-3 

1 

-23.5 

-10.25 

62 

24.5 

0 

2 

-21 

-10.25 

63 

24.5 

3 

3 

-18 

-10.25 

64 

24.5 

6 

4 

-15 

-10.25 

65 

24.5 

9 

5 

-12 

-10.25 

66 

24.5 

12.25 

6 

-9 

-10.25 

Inside  of  upper  insulation  panel 

7 

-6 

-10.25 

67 

22.5 

13 

8 

-3 

-10.25 

68 

21 

13 

9 

0 

-10.25 

69 

18 

13 

10 

3 

-10.25 

70 

15 

13 

11 

6 

-10.25 

71 

12 

13 

12 

9 

-10.25 

72 

9 

13 

13 

12 

-10.25 

73 

6 

13 

14 

15 

-10.25 

74 

3 

13 

15 

18 

-10.25 

75 

0 

13 

16 

21 

-10.25 

76 

-3 

13 

17 

23.5 

-10.25 

77 

-6 

13 

Outside  of  bottom  insulation  panel 

78 

-9 

13 

18 

-23.5 

-10.75 

79 

-12 

13 

19 

-21 

-10.75 

80 

-15 

13 

20 

-18 

-10.75 

81 

-18 

13 

21 

-15 

-10.75 

82 

-21 

13 

22 

-12 

-10.75 

83 

-22.5 

13 

23 

-9 

-10.75 

Outside  of  upper  insulation  panel 

24 

-6 

-10.75 

84 

22.5 

13.5 

25 

-3 

-10.75 

85 

21 

13.5 

26 

0 

-10.75 

86 

18 

13.5 

27 

3 

-10.75 

87 

15 

13.5 

28 

6 

-10.75 

88 

12 

13.5 

29 

9 

-10.75 

89 

9 

13.5 

30 

12 

-10.75 

90 

6 

13.5 

31 

15 

-10.75 

91 

3 

13.5 

32 

18 

-10.75 

92 

0 

13.5 

33 

21 

-10.75 

93 

-3 

13.5 

34 

23.5 

-10.75 

94 

-6 

13.5 

Inside  of  left  insulation  panel 

95 

-9 

13.5 

35 

-24 

-9.25 

96 

-12 

13.5 

36 

-24 

-6 

97 

-15 

13.5 

37 

-24 

-3 

98 

-18 

13.5 

38 

-24 

0 

99 

-21 

13.5 

39 

-24 

3 

100 

-22.5 

13.5 

40 

-24 

6 

101 

12  o'clock  on  right  pipe  (8  in.) 

41 

-24 

9 

102 

3  o'clock  on  right  pipe  (8  in.) 

42 

-24 

12.25 

103 

6  o'clock  on  right  pipe  (8  in.) 

Outside  of  left  insulation  panel 

104 

9  o'clock,  right  pipe  (8  in.) 

43 

-24.5 

-9.25 

105 

12  o'clock,  right  pipe  insulation 

44 

-24.5 

-6 

106 

3  o'clock,  right  pipe  insulation 

45 

-24.5 

-3 

107 

6  o'clock,  right  pipe  insulation 

46 

-24.5 

0 

108 

9  o'clock,  right  pipe  insulation 

47 

-24.5 

3 

109 

12  o'clock  on  left  pipe  (4  in.) 

48 

-24.5 

6 

110 

9  o'clock  on  left  pipe  (4  in.) 

49 

-24.5 

9 

111 

6  o'clock  on  left  pipe  (4  in.) 

50 

-24.5 

12.25 

112 

3  o'clock  on  left  pipe  (4  in.) 

Inside  of  right  insulation  panel 

113 

12  o'clock,  left  pipe  insulation 

51 

24 

-9.5 

114 

9  o'clock,  left  pipe  insulation 

52 

24 

-6 

115 

6  o'clock,  left  pipe  insulation 

53 

24 

-3 

116 

3  o'clock,  left  pipe  insulation 

54 

24 

0 

In  air,  89  in. 

from  end 

55 

24 

3 

117 

-18 

-1 

56 

24 

6 

118 

-6 

-1 

57 

24 

9 

119 

6.5 

-1 

58 

24 

12.25 

120 

20 

-1 

Outside  of  right  insulation  panel 

121 

coolant  supply 

59 

24.5 

-9.5 

122 

coolant  return 

60 

24.5 

-6 

29 
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Figure  20.  Data  acquisition  system  for  use  with  the  1-ft  x  1-ft  enclosure. 


and  y  coordinates;  these  were  placed  66  in.  from  the  coolant  supply  end,  except 
where  noted.  The  total  interior  length  was  106  in. 

Procedure 

Once  the  apparatus  was  connected  to  a  coolant  supply  and  electrical  power 
source,  several  adjustments  were  possible.  The  coolant  temperature  could  be 
adjusted,  the  coolant  flow  through  each  panel  could  be  controlled,  and  the  energy 
input  to  the  pipe  heater(s)  could  also  be  adjusted.  Experiments  consisted  of  obtain¬ 
ing  a  range  of  temperature  conditions  for  both  the  interior  walls  of  the  apparatus 
and  pipe(s)  at  steady-state  conditions.  Steady  state  was  determined  to  be  obtained 
once  the  temperatures  along  the  top  of  the  enclosure  were  changing  by  a  small 
amount  and  in  an  apparently  random  fashion.  Once  a  steady-state  condition  was 
reached,  three  sets  of  data  from  all  of  the  thermocouples  and  the  watt-meter  were 
stored  on  an  approximately  hourly  basis.  Thus  for  each  condition,  three  sets  of  data 
were  collected.  Temperature  data  from  the  2-ft  x  4-ft  enclosure  were  collected  con¬ 
tinually  on  an  hourly  basis. 


RESULTS 

General  information 

Data  obtained  from  the  experimental  and  numerical  experiments  are  summa¬ 
rized  in  Appendixes  A  and  B,  respectively.  Physical  data  are  based  on  the  average 
of  three  hourly  readings.  Temperatures,  Nusselt  number,  Rayleigh  number,  heat 
flow  though  the  interior  wall,  and  the  average  thermal  conductance  at  the  interior 
enclosure  surface  are  presented.  Heat  flow  through  each  side,  assuming  1  ft  of 
enclosure  length,  was  calculated  by  using  the  temperatures  around  the  insulation 
and  averaging  the  heat  flows  calculated  at  each  thermocouple  location  (physical 
experiments)  or  the  temperatures  at  each  node  location  (numerical).  These  values 
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were  averaged  to  obtain  a  heat  flow  value  for  each  side  (top,  bottom,  left,  and  right) 
and  an  overall  average  value.  Nusselt  and  Rayleigh  numbers  were  then  calculated 
using  the  surface  conductances  obtained  using  the  average  interior  surface  tem¬ 
perature  for  material  properties  and  the  temperature  difference  between  the  aver¬ 
age  pipe  and  average  insulation  surface  temperatures,  for  each  side  and  for  the 
overall  average. 

The  following  equations  developed  from  data  in  Raznjevic  (1976)  were  used  for 
calculating  material  properties. 

Viscosity  of  air  (ft/s2) 

D  =  1.27573E  -  04+6.1411E-  07TAVG.  (126) 

Density  of  air  (lbm/ft3) 


pair  =  8.42416E  -  02  - 1 .93863E  -  04TAVG  +  4.16195E  -  07TjVG  •  (127) 

Volumetric  specific  heat  of  air  (Btu/ft3oF) 

Cv=0.241pair.  (128) 

Coefficient  of  thermal  expansion  (1/°F) 

P  =  2.177E-  03  -  4.74865E-  06  TAVG  +  9.42743E-  09T^vg 
-1.04328E-11T*vg  (129) 

Thermal  conductivities  (Btu/fthr°F) 

fcair  =  0.01309  +  2.14766E  -  05TAVG  (130) 

fcEPS  =  (0.214583  +  5.36829E  -  04TAVG )  / 12  (131) 

fcpipe  ^ul  =  (O.196851e°-00211687TAVG )  / 12 .  (132) 


TAvg  is  the  average  temperature  of  the  two  pertinent  surfaces,  i.e.,  for  air  the 
average  is  that  of  the  inside  EPS  surface  and  the  pipe  or  pipe  insulation  surface 
temperatures. 

The  hypothetical  gap  width  was  used  as  the  length  parameter  in  the  Nusselt 
and  Rayleigh  number  calculations;  these  values  and  other  enclosure  dimensions 
are  in  Table  8.  Table  9  shows  the  physical  configurations  and  radiation  emissivity 
values  for  the  numerical  experiments.  Three  sets  of  emissivity  values  were  used; 
two  different  values  for  EPS  were  chosen  to  represent  new  (0.6)  and  old  (0.9)  insu¬ 
lation.*  For  the  pipe  and  pipe  insulation  the  two  values  were  for  aluminized  paint 
(0.5)  and  no  paint  (0.9). 


*  Personal  communication  ,  Stephen  N.  Flanders,  U.S.  Army  Cold  Regions  Research  and  Engineer¬ 
ing  Laboratory,  1994. 
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Table  8.  Enclosure  dimensions  and  effective  gap  widths. 


Outside  radius  of  Effective  radius 


Configuration 

Pipe  description 

pipe/insulation 

of  enclosure 

Effective  gap 

1 

4-in.  bare 

0.18750  ft 

0.583568  fta 

0.396068  ft 

2 

4-in.  insulated 

0.27083 

0.583568 

0.312738 

3 

2-in.  bare 

0.09896 

0.583568 

0.484608 

4 

2-in.  insulated 

0.18229 

0.583568 

0.401278 

5 

4-in.  insulatedb 

0.35416 

0.755456c 

0.401296 

6 

2-in.  insulated 

0.18229 

1.220187d 

1.037897 

7 

2-  and  2-in.  insulated 

0.364583 

1.220187d 

0.855605 

8 

4-  and  8-in.  insulated 

0.713542e 

1.856808f 

1.143266e 

a  1-ft  x  1-ft  enclosure 
b  2  in.  of  insulation 
c  1.27-ft  x  1.27-ft  enclosure 
d  2-ft  x  2-ft  enclosure 

e  0.442708  and  1.4140,  if  only  8-in.  pipe  heated 
f  2-ft  x  4-ft  enclosure 


Table  9.  Configurations  for  the  numerical  experiments. 


Enclosure  outside 
dimensions 

Number  of 
pipes 

Pipe 

diameter(s) 

Pipe  insulation 
thickness 

Emissivity  values 
Pipe 1  Insulation 

1  ft  X  1  ft 

1 

4.5  in. 

0.0  in. 

0.9 

0.9 

1 

4.5 

0.0 

0.9 

0.6 

1 

4.5 

0.0 

0.5 

0.9 

1 

4.5 

0.0 

No  radiation 

1 

4.5 

1.0 

0.9 

0.9 

1 

4.5 

1.0 

0.9 

0.6 

1 

4.5 

1.0 

0.5 

0.9 

1 

2.375 

0.0 

0.9 

0.9 

1 

2.375 

0.0 

0.9 

0.6 

1 

2.375 

0.0 

0.5 

0.9 

1 

2.375 

1.0 

0.9 

0.9 

1 

2.375 

1.0 

0.9 

0.6 

1 

2.375 

1.0 

0.5 

0.9 

1 

2.375 

1.0 

No  radiation 

1.27  ft  x  1.27  ft 

1 

4.5 

2.0 

0.9 

0.9 

2  ft  x  2  ft 

1 

2.375 

1.0 

0.9 

0.9 

2 

2.375, 2.375 

1.0, 1.0 

0.9 

0.9 

2  ft  x  4  ft 

2 

4.5,  8.625 

1.0, 1.0 

0.9 

0.9 

1  or  pipe  insulation 


Experimental  data 

The  complete  set  of  temperature  and  power  measurements  is  reported  in  Rich¬ 
mond  et  al.  (1997).  The  average  Nusselt  and  Rayleigh  numbers  obtained  for  the 
physical  experiments  conducted  with  the  1-ft  x  1-ft  enclosure  are  plotted  in  Fig¬ 
ures  21-24.  Vertical  scatter  is  due  to  the  radiation  effect  on  the  surface  conductance; 
a  similar  effect  is  seen  in  the  numerical  data,  shown  in  the  next  section.  In  Figure 
21,  the  1991  data  result  in  slightly  higher  Nusselt  numbers.  Some  of  this  difference 
is  due  to  the  temperature  effect  on  the  conductance,  but  some  may  also  be  due  to 
the  larger  thermocouples  used  to  measure  surface  temperatures.  The  effect  of  paint¬ 
ing  the  pipe  with  the  aluminum  paint,  with  a  subsequent  lower  emissivity,  is  clearly 
seen,  resulting  in  lower  Nusselt  numbers  compared  to  the  unpainted  pipe.  Figure 
22  shows  similar  scatter  with  apparently  one  significant  outlier,  which  occurred  in 
the  data  obtained  on  28  September  1996  (see  App.  A).  By  plotting  the  calculated 
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Figure  21.  Nusselt  and  Rayleigh  number  plot  for  the  4-in.  pipe  in  the 
1-ft  xl-ft  enclosure  (experimental  data). 


105  io6 

Rayleigh  Number 


Figure  22.  Nusselt  and  Rayleigh  number  plot  for  the  insulated  4-in. 
pipe  in  the  1-ft  xl-ft  enclosure  (experimental  data). 

heat  flux  data  versus  the  difference  between  the  average  interior  wall  temperature 
and  the  average  pipe  temperature,  it  was  found  that  a  number  of  data  points  did 
not  agree  with  a  general  linear  trend  of  the  data.  This  was  eventually  traced  to  events 
involving  the  low  temperature  coolant  supply  and  the  fact  that  the  low  tempera¬ 
ture  chiller  was  being  run  manually  In  these  cases  the  coolant  became  too  warm 
before  an  operator  was  able  to  get  the  chiller  back  on  line.  This  was  not  observed  in 
the  thermocouple  data  being  monitored;  although  it  appeared  that  steady  condi¬ 
tions  had  been  reached,  apparently  temperatures  were  still  changing  slowly 
Figure  23  contains  the  data  for  the  2-in.  uninsulated  pipe;  two  outliers  are 
observed.  Figure  24  contains  the  data  for  the  2-in.  insulated  pipe;  two  outliers  are 
in  this  plot  although  they  cannot  easily  be  discerned.  Those  data  determined  to  be 
outliers  are  indicated  in  Appendix  A  and  were  not  used  in  further  analysis  or  com¬ 
parisons. 
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Figure  23.  Nusselt  and  Rayleigh  number  plot  for  the  2-in.  pipe  in  the 
1-ft  xl-ft  enclosure  (experimental  data). 
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Figure  24.  Nusselt  and  Rayleigh  number  plot  for  the  2-in.  insulated  pipe 
in  the  1-ft  xl-ft  enclosure  (experimental  data). 


Experimental  data  obtained  for  the  2-ft  x  4-ft  enclosure  are  somewhat  limited. 
Time  constraints  resulted  in  only  one  pipe  configuration  and  limited  temperature 
combinations.  Steady-state  temperatures  took  longer  to  achieve  (compared  to  the 
1-ft  x  1-ft  enclosure),  and  even  then  may  have  been  influenced  by  the  room  tem¬ 
perature,  which  fluctuated  over  a  ±5°F  temperature  range.  Some  low  temperature 
tests  were  attempted  with  only  the  8-in.  pipe  heated;  however,  the  pipe  heater  was 
not  able  to  hold  the  desired  temperature  (at  or  above  235°F)  for  all  desired  tests. 
Figure  25  shows  the  Nusselt  and  Rayleigh  number  plots.  Gap  width  was  deter¬ 
mined  by  using  an  effective  radius  of  the  heated  pipe(s);  the  values  are  in  Table  8. 
The  interior  temperature  range  for  these  data  is  from  about  3°F  to  80°F.  This 
almost  spans  the  range  of  interest  for  most  utilidor  designs,  even  though  a  rather 
narrow  range  of  Rayleigh  numbers  was  obtained. 
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Figure  25.  Nusselt  and  Rayleigh  number  plot  for  the  2  -ft  x4-ft  enclosure 
(experimental  data). 

There  is  significant  difference  in  the  curve  shapes  between  the  two-pipe  heat¬ 
ing  conditions.  This  may  be  due  to  more  stratification  of  the  air  in  the  bottom  of 
the  enclosure  combined  with  a  larger  temperature  difference  compared  to  the  single 
heated  pipe  configuration;  this  affects  the  heat  conductance  (h)  in  the  Nusselt  num¬ 
ber.  Comparing  Nusselt  numbers  for  the  bottom  surfaces  (App.  A),  it  can  be  seen 
that  there  is  a  greater  change  for  two  heated  pipes  compared  to  the  single  heated 
pipe.  This  observation  also  explains  the  curved  shape  of  the  Nu~Ra  number  data 
obtained  with  numerical  data  discussed  in  the  next  section. 

Numerical  data 

The  finite-element  computer  program  FECOME,  described  earlier,  was  used  to 
obtain  additional  heat  transfer  data  from  numerical  experiments.  The  objective  of 
the  numerical  experiments  was  to  extend  the  database  of  enclosure  configurations 
and  boundary  conditions,  and  to  make  comparisons  with  the  physical  experiments. 
The  numerical  experiments  allowed  calculations  to  be  made  without  radiation 
boundaries  and  with  different  combinations  of  emissivity  values. 

Figure  26  shows  one  of  the  meshes  used 
for  the  uninsulated  4-in.  pipe  in  the  1-ft  x  1- 
ft  enclosure.  It  has  the  same  internal  dimen¬ 
sions  as  the  experimental  apparatus,  includ¬ 
ing  the  0.5-in.  layer  of  EPS  insulation 
enclosing  the  cavity  of  air.  Temperatures 
around  the  outside  of  the  insulation  were 
held  constant  and  the  surface  representing 
the  outer  diameter  of  the  pipe  was  held  at  a 
series  of  temperatures.  Increasing  Rayleigh 
number  values  were  obtained  by  decreasing 
the  outside  boundary  temperatures  in  5-  or 
10-degree  increments  until  FECOME  was  no 
longer  able  to  converge  to  a  solution.  A  small 
temperature  change  in  boundary  condi¬ 
tions,  and  the  use  of  a  previous  solution  as 


Figure  26.  Mesh  for  the  1-ft  xl-ft 
enclosure  with  a  4-in.  pipe  (3,552 
nodes ,  1,152  elements,  10,848  d.o.f). 
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Figure  27.  Mesh  for  the  2  -ft  x2-ft  enclo¬ 
sure  with  a  2-in.  insulated  pipe  (16,704 
nodes,  5,972  elements,  48,862  d.o.f ). 
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an  initial  solution  estimate,  aided  the 
model  in  converging  to  a  solution. 

Towards  the  end  of  this  study  it  was  found 
that  often  FECOME  converged  to  an 
oscillatory  solution,  which  may  actually 
occur  in  steady  solutions.  The  tempera¬ 
tures  were  generally  within  or  close  to  the 
convergence  criteria,  while  the  maximum 
velocity  and  pressure  changes  were  small 
(2-3%).  In  Appendix  B,  the  solutions  that 
were  oscillating  are  indicated  with  an 
asterisk  in  the  file  name. 

The  effect  of  mesh  density  was  investi¬ 
gated  and  reported  in  Richmond  (1997).  In 
general,  it  was  found  that  mesh  density 
had  little  effect  on  average  values,  but 
denser  meshes  were  required  to  obtain 

solutions  at  higher  Rayleigh  numbers .  The  meshes  generated  for  all  of  the  1  -ft  square 
enclosure  configurations  are  similar  to  Figure  26,  as  is  the  mesh  for  the  1.27-ft  x 
1.27-ft  enclosure.  Meshes  for  the  other  configurations  are  shown  below  (Fig.  27- 

29).  A  limited  number  of 
solutions  were  obtained  for 
the  2-ft  x  2-ft  enclosure  with 
two  pipes  but  none  for  the 
2-ft  x  4-ft  enclosure.  The  rea¬ 
son  for  this  is  not  clear,  as 
these  meshes  are  as  dense  as 
those  used  for  the  1-ft  x  1-ft 
enclosures.  Attempts  to 
solve  these  configurations 
required  large  quantities  of 
memory  and  cpu  time.  It 
appears  that  the  solution 
methods  used  by  FECOME  may  not  be 
optimal  for  these  large-scale  problems. 
Occasionally  reports  have  been  made  in 
the  literature  regarding  difficulties  in 
obtaining  solutions  to  high  (107)  Ray¬ 
leigh  number  problems,  but  there  has 
been  no  indication  of  the  best  way  to 
solve  this  problem. 

Figures  30-33  display  the  numerical 
data  in  Ra-Nu  number  format;  dashed 
lines  connect  data  obtained  with  the 
same  pipe  temperatures  (values  of  250, 
150,  100,  80,  and  40°F).  Data  are  also 
shown  for  the  conditions  where  no 
radiation  was  modeled.  The  three  dif¬ 
ferent  combinations  of  radiation  emis- 
sivities  are  designated  as  follows: 


Figure  28.  Mesh  for  the  2-ft  x4-ft  enclosure , 
8-in.  insulated  pipes  (13,338  nodes ,  4,356 
38,834  d.o.f). 
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Figure  29.  Mesh  for  the  2-ft  x2-ft  enclosure 
with  two  2-in.  pipes  (14,759  nodes,  4,824 
elements,  42,188  d.o.f). 
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Figure  30.  Nusselt  and  Rayleigh  number  plot  for  the  4-in.  pipe  in  the 
1-ft  xl-ft  enclosure  (numerical  data). 


pipe  in  the  1-ft  xl-ft  enclosure  (numerical  data). 

Radiation  A,  Pipe  or  insulation  surface  -0.9,  inside  insulation  surface  -0.9 
Radiation  B,  Pipe  or  insulation  surface  -0.9,  inside  insulation  surface  -0.6 
Radiation  C,  Pipe  or  insulation  surface  -0.5,  inside  insulation  surface  -0.9. 

Significant  differences  are  seen  for  each  different  radiation  condition. 

In  Figure  30,  the  convection  data  alone  (no  radiation  condition)  do  not  all  coin¬ 
cide  on  one  line  as  expected.  This  is  thought  to  be  due  to  errors  associated  with 
too  coarse  a  mesh,  or  for  other  unidentified  reasons. 

Figure  34  shows  the  data  for  the  1.27-ft  x  1.27-ft  and  the  2-ft  x  2-ft  enclosures. 
The  1.27-ft  x  1.27-ft  enclosure  had  a  4-in.  pipe  with  two  inches  of  insulation.  This 
resulted  in  nearly  the  same  effective  gap  as  the  2-in.  insulated  pipe  in  the  1-ft  x  1- 
ft  enclosure.  The  meager  amount  of  data  obtained  with  two  pipes  is  also  shown 
on  the  plot.  The  effective  gap  in  this  case  was  obtained  using  the  total  area  of  the 
two  insulated  pipes  to  determine  an  effective  interior  radius. 
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Figure  32 .  Nusselt  and  Rayleigh  number  plot  for  the  2-in.  pipe  in 
the  1-ft  xl-ft  enclosure  (numerical  data). 


Figure  33.  Nusselt  and  Rayleigh  number  plot  for  the  insulated  2-in. 
pipe  in  the  1-ft  xl-ft  enclosure  (numerical  data). 


Figure  34.  Nusselt  and  Rayleigh  number  plot  for  the  1.27 -ft  xl.27 - 
ft  and  2 -ft  x2-ft  enclosures  (numerical  data). 


Go  to  Contents 


38 


Figure  35.  Heat  flux  from  the  4-in.  pipe  through  the  1-ft  xl-ft  enclosure. 


ANALYSIS 

The  Nusselt  Number-Rayleigh  Number  plots  in  the  previous  section  showed 
that  no  simple  direct  correlation  between  these  parameters  would  be  found.  Com¬ 
parisons  of  other  parameters  are  made  in  this  section  and  a  new  approximation 
for  the  effective  conductivity  of  air  is  proposed.  Comparisons  of  numerical  solu¬ 
tions  using  this  effective  conductivity  correlation  are  made  with  those  obtained 
using  FECOME  and  with  experimental  data  from  the  2-ft  x  4-ft  enclosure. 

Figures  35-38  compare  heat-flux-per-foot  data  from  the  numerical  and  experi¬ 
mental  results  with  the  eight  methods  presented  in  Table  1  for  four  configurations 
of  the  1-ft  x  1-ft  enclosure.  Figures  35  (a)  and  (b)  compare  the  effect  of  emissivity 
values  for  radiation  conditions  A  and  C.  In  these  figures,  temperature  difference 
is  the  total  temperature  difference  for  the  system;  heat  flux  is  through  the  mean 
perimeter  of  the  enclosure.  Using  the  methods  in  Table  1  to  determine  the  effec¬ 
tive  conductivity  of  the  air,  the  heat  flux  was  calculated  using  eq  51-57.  Vertical 
scatter  within  a  given  method  is  related  to  the  average  interior  temperatures.  The 
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Figure  36.  Heat  flux  from  the  2-in.  pipe  through  the  1-ft  xl-ft  enclosure. 


Temperature  Difference  (°F) 


Figure  37.  Heat  flux  from  the  4-in.  insulated  pipe  through  the  1-ft  xl- 
ft  enclosure. 
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Figure  38.  Heat  flux  from  the  2-in.  insulated  pipe  through  the  1-ft  xl- 
ft  enclosure. 
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numerical  and  experimental  data  compare  well  in  all  five  plots*  For  the  uninsulated 
pipes,  the  best  agreement  between  the  numerical  and  experimental  data  is  with 
methods  1  and  2,  which  use  eq  37  and  38  to  calculate  the  effective  conductivity. 
The  differences  caused  by  emissivity  values  can  also  be  seen,  with  a  slightly 
reduced  heat  flux  observed  for  the  lower  emissivity  (condition  C).  For  the  insulated 
pipes,  there  is  general  agreement  between  all  the  methods.  This  occurs  because  the 
thermal  resistance  due  to  the  air  gap  becomes  small  relative  to  the  resistance  of  the 
enclosure  and  pipe  insulation. 

Figures  39-42  compare  the  ratio  of  effective  conductivity  to  the  thermal  conduc¬ 
tivity  of  air  (fceff/ta ir)  with  the  average  temperature  of  the  interior  surfaces.  The 
effective  conductivity  was  calculated  using  eq  36  and  the  heat  flow  through  the 
inside  surface.  Good  agreement  is  seen  between  the  numerical  and  experimental 
test  data,  while  comparison  of  the  Table  1  methods,  in  some  cases,  shows  effects  of 
temperature  difference.  Some  of  the  methods  show  the  same  trend  as  the  numerical 
and  experimental  data  but  differ  in  magnitude  (for  example,  method  4  in  Fig.  42). 


Figure  39.  Ratio  of  effective  conductivity  to  the  conductivity  of  air  vs. 
the  average  interior  temperature  (4-in.  pipe). 
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Figure  40.  Ratio  of  effective  conductivity  to  the  conductivity  of  air 
vs.  the  average  interior  temperature  (2-in.  pipe). 
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Figure  41.  Ratio  of  effective  conductivity  to  the  conductivity  of  air  vs. 
the  average  interior  temperature  (4-in.  insulated  pipe). 
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Figure  42.  Ratio  of  effective  conductivity  to  the  conductivity  of  air  vs. 
the  average  interior  temperature  (2-in.  insulated  pipe). 


Curves  of  the  following  form  were  fit  to  the  numerical  and  experimental  data 
in  Figures  39-42,  and  to  the  additional  numerical  data: 

=  AeBTAVG  (133) 

k  . 

*air 

where  TAVG  is  the  average  of  the  interior  surface  temperatures  and  A  and  B  are 
defined  in  Table  10.  Equations  134-140  are  plotted  in  Figure  43.  Comparing  eq  134 
with  eq  139  and  140  shows  a  reduction  of  39%  and  15%  from  radiation  conditions 
A  to  B  and  A  to  C,  respectively. 

In  an  attempt  to  correlate  the  coefficients  in  eq  134-138  with  a  geometric  param¬ 
eter  associated  with  the  enclosure,  it  was  found  that  a  slight  linear  correlation 
exists  between  the  radius  (r)  of  the  interior  pipe  (or  insulation)  and  the  parameter 
A  (the  intercept).  These  data  and  the  correlation  are  shown  in  Figure  44.  Using  an 
average  value  of  the  slopes  (B)  results  in  the  following  equation: 
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Table  10.  Coefficients  for  eq  133. 


Pipe  or 
insulation 

Effective  gap  radius  A  B  eq  Description 


0.396068 

0.1875 

12.026 

0.003094 

134 

0.312738 

0.27083 

10.9327 

0.003057 

135 

0.484608 

0.9896 

9.7324 

0.004123 

136 

0.401278 

0.18229 

12.4199 

0.001706 

137 

0.401296 

0.35416 

13.2964 

0.003680 

138 

0.396068 

0.1875 

7.205 

0.003387 

139 

0.396068 

0.1875 

10.059 

0.003424 

140 

4-in.  pipe  in  the  1-ft  x  1-ft  enclosure,  numerical  and 
experimental  data,  emissivities:  pipe,  0.9;  enclosure,  0.9. 
4-in.  insulated  pipe  in  the  1-ft  x  1-ft  enclosure,  numeri¬ 
cal  and  experimental  data  emissivities:  pipe  insulation, 
0.9;  enclosure,  0.9. 

2-in.  pipe  in  the  1-ft  x  1-ft  enclosure,  numerical  and 
experimental  data,  emissivities:  pipe,  0.9;  enclosure,  0.9. 
2-in.  insulated  pipe  in  the  1-ft  x  1-ft  enclosure,  numeri¬ 
cal  and  experimental  data,  emissivities:  pipe  insula¬ 
tion,  0.9;  enclosure,  0.9. 

4-in.  pipe  with  2  in.  of  insulation  in  the  1.27-ft  x  1.27-ft 
enclosure,  emissivities:  pipe  insulation,  0.9;  enclo¬ 
sure,  0.9. 

4-in.  pipe  in  the  1-ft  x  1-ft  enclosure,  numerical  and 
experimental  data,  emissivities:  pipe,  0.9;  enclosure,  0.6. 
4-in.  pipe  in  the  1-ft  x  1-ft  enclosure,  numerical  and 
experimental  data,  emissivities:  pipe,  0.5;  enclosure,  0.9. 


Figure  43.  Effective  conductivity  correlations. 


Figure  44.  Correlation  of  pipe  radius  with  intercepts. 
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(141) 


^  =  (9.5031+  9.9585r)ea00373TAVG . 

*air 

Table  11  presents  comparisons  at  three  air  temperatures  for  eq  134  through  138, 
and  141.  The  poorest  comparisons  are  with  eq  135  and  137,  which  have  differences 
of  -19  and  17%  at  an  air  temperature  of  100°F. 

The  results  from  three  numerical  experiments  were  next  compared  with  corre¬ 
sponding  numerical  experiments  in  which  the  conductivity  of  air  was  specified 
using  eq  141  and  treated  as  a  solid  (without  radiation);  the  same  FE  meshes  were 
used  in  both  cases.  Table  12  contains  the  parameters  and  description  of  each  com¬ 
parison.  Figures  45-50  compare  the  inside  surface  temperatures  and  temperature 
contours  for  each  comparison. 

In  the  figures  above,  it  can  be  seen  that  using  an  effective  conductivity  in  lieu  of 
the  full  convection  solution  will  produce  inaccurate  temperature  distributions  in 
the  air.  In  some  cases,  the  predicted  interior  surface  temperatures  agree,  but  the 
quality  of  agreement  depends  on  which  surface  (top,  bottom,  or  side)  is  being  con¬ 
sidered.  Better  agreement  is  seen  when  a  small  temperature  difference  exists  within 
the  "air,"  and  when  the  pipe  is  insulated.  Average  inside  surface  temperatures 


Table  11.  Comparison  of  effective  conductivity 
correlations. 


Pipe  radius  (ft) 


Iavg 

0.1875 

0.27083  0.09896 

0.18229 

0.35416 

Equation  141,  k{ 

eff^air 

20 

12.25 

13.15 

11.30 

12.20 

14.04 

60 

14.22 

15.26 

13.12 

14.16 

16.30 

100 

16.51 

17.72 

15.23 

16.44 

18.92 

Equations,  fceff//cair 

134 

135 

136 

137 

138 

20 

12.79 

11.62 

10.57 

13.65 

14.31 

60 

14.48 

13.13 

12.46 

16,47 

16.58 

100 

16.39 

14.84 

14.70 

19.88 

19.21 

Residuals 

20 

0.543 

-1.523 

-0.732 

1.450 

0.273 

60 

0.257 

-2.127 

-0.655 

2.314 

0.284 

100 

-0.124 

-2.874 

-0.531 

3.448 

0.291 

%  differences 

20 

4.2 

-13.1 

-6.9 

10.6 

1.9 

60 

1.8 

-16.2 

-5.3 

14.1 

1.7 

100 

-0.8 

-19.4 

-3.6 

17.3 

1.5 

Table  12.  Conditions  for  comparisons  between  numerical  conduc¬ 
tion  and  convection  solutions. 


Name* 

Pipe 

temp.  ( °F) 

Outside 
temp.  ( <P) 

Pipe  radius  (ft) 

Estimated 

T Avg<*> 

K/f 

( Btu/hrft 2  7) 

sq4ib65a 

150 

85 

0.27083 

100 

0.26995 

sq2d35c 

80 

45 

0.09896 

60 

0.18865 

sq2ic65a 

150 

85 

0.18229 

100 

0.25045 

*Names  correspond  to  file  names  in  Appendix  B. 
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Figure  45.  Temperature  contours  for 
sq4ib65a.  Solid  lines  are  from  the  convec¬ 
tion  and  radiation  solution ,  dashed  lines  are 
from  the  conduction  solution. 


Distance  Around  Inside  Insulation  Surface 
(Degrees  from  Vertical) 

Figure  46.  Inside  surface  temperatures  for  sq4ib65a. 


Figure  47.  Temperature  contours  for 
sq2d35c.  Solid  lines  are  from  the  convection 
and  radiation  solution ,  dashed  lines  are 
from  the  conduction  solution. 
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Distance  Around  Inside  Insulation  Surface 
(Degrees  from  Vertical) 

Figure  48.  Inside  surface  temperatures  for  sq2d35c. 


Figure  49.  Temperature  contours  for 
sq2ic65a.  Solid  lines  are  from  the  convec¬ 
tion  and  radiation  solution ,  dashed  lines 
are  from  the  conduction  solution. 


Distance  Around  Inside  Insulation  Surface 
(Degrees  from  Vertical) 

Figure  50.  Inside  surface  temperatures  for  sq2ic65a. 
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Figure  51.  Ratio  of  effective  conductivity  to  the  conductivity  of  air  vs. 
the  average  interior  temperature  (2-ft  x4-ft  enclosure). 


for  the  three  cases  were  as  follows:  sq4ib65a,  95.20°F  and  95.20°F;  sq2d35c,  53.56°F 
and  55.25°F;  sq2ic65a,  91.83°F  and  90.82°F,  for  the  convection  versus  conduction 
solutions,  respectively.  These  average  values  agree  very  well,  and  if  used  in  calcu¬ 
lations  of  average  heat  loss,  would  give  comparable  values.  Examining  the  insu¬ 
lation  enclosure  temperature  contours,  it  can  be  seen  that  approximately  midway 
through  the  insulation  the  temperature  contours  agree  fairly  well,  but  one  reason 
for  this  agreement  is  the  fixed  outer  boundary  temperatures. 

A  similar  approach  was  followed  using  the  2-ft  x  4-ft  experimental  data.  Figure 
51  shows  the  comparison  between  the  methods  from  Table  1  and  the  experimen¬ 
tal  data  obtained  with  both  pipes  heated;  weighted  averages  of  the  pipe  insula¬ 
tion  surface  temperatures  were  used  (note  that  the  Table  1  methods  are  for  two 
pipes  heated  at  the  same  temperature,  and  shouldn't  be  compared  with  the 
experimental  data  for  the  single  heated  pipe).  These  values  are  much  higher  than 
those  obtained  from  the  smaller  enclosure,  and  the  intercepts  do  not  correlate  with 
pipe  (insulation)  radius  as  determined  earlier.  The  slopes  are  near  zero,  and  the 
intercepts  (from  curve  fits)  are  56.5836  and  86.6341  for  the  cases  of  both  pipes  heated 
and  a  single  pipe  heated,  respectively. 

Because  no  numerical  convection  and  radiation  solutions  were  obtained  for  this 
size  enclosure,  comparisons  could  be  made  only  with  the  experimental  data. 
Using  the  physical  experimental  data  from  13  January  1997  and  FECOME  with  an  ef¬ 
fective  conductivity  determined  by  two  methods,  comparisons  were  made  between 
predicted  and  measured  interior  enclosure  surface  temperatures.  Polynomial  curve 
fits  were  made  to  the  temperatures  measured  on  the  exterior  of  the  0.5-in.  insula¬ 
tion  (Fig.  52)  and  were  applied  as  fixed  boundary  temperatures  to  the  exterior  sur¬ 
face  nodes  of  the  mesh  shown  in  Figure  53.  The  average  pipe  temperatures,  142.75 °F 
and  146.20°F,  were  used  for  the  4-  and  8-in.  pipes,  respectively.  Two  values  for  the 
effective  conductivity  of  air  were  used.  One  value  was  determined  by  using  the 
effective  radius  of  the  two  pipes  (0.71354  ft)  and  eq  141,  which  yields  a  fceff  of 
0.22863.  The  other  value  was  0.752115,  determined  from  fceff//cair = 56.584  (the  curve 
fit  to  the  test  data  in  Fig.  51),  which  resulted  in  a  ke #  of  0.1773.  The  air  temperature 
(9.39°F)  from  the  experimental  test  data  was  used  to  determine  fcair. 
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Location  on  Given  Surface  (ft) 

Figure  52.  Polynomial  curve  fits  to  experimental  boundary  data. 


Figure  53.  Finite  element  mesh  used  for  conduction  solution 
of  the  2-ft  x4-ft  enclosure. 


Figure  54.  Comparison  of  experimental  and  numerical  inside  insula¬ 
tion  surface  temperatures  (data  from  13  January  1997). 
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Figure  54  compares  the  measured  and  calculated  interior  insulation  surface 
temperatures.  Fair  agreement  is  obtained  for  all  surfaces  except  the  bottom, 
where  temperatures  are  much  warmer  than  measured.  It  can  also  be  observed  that 
the  temperatures  are  not  very  sensitive  to  changes  in  the  effective  conductivity. 
The  case  chosen  was  that  which  had  the  greatest  temperature  difference  across  the 
air  gap,  and  this  should  be  considered  when  comparing  the  temperatures.  It  could 
also  be  noted  that  closeness  of  the  outside  boundary  conditions  to  the  compared 
temperatures  are  ensuring  reasonable  results.  In  defense  of  this,  in  an  actual  utilidor 
design  the  wall  thickness  would  be  much  greater,  but  with  well-known  thermal 
properties.  Because  conduction  solutions  can  be  obtained  with  a  high  degree  of 
confidence,  then  temperatures  in  similar  locations  will  in  most  cases  be  known 
fairly  well. 


SUMMARY 

Three  approaches  to  the  thermal  analysis  of  utilidors  were  investigated:  the  tra¬ 
ditional  or  currently  accepted  practice  of  one-dimensional  analysis,  numerical 
analysis  with  modeling  of  convection  and  radiation,  and  numerical  conduction 
analysis  using  an  effective  conductivity  to  account  for  convection  and  radiation 
effects.  Each  method  has  limitations  and  advantages. 

The  one-dimensional  analysis  did  not  produce  good  results  when  uninsulated 
pipes  were  modeled;  only  some  correlations  can  account  for  multiple  pipes,  off- 
center  locations,  or  other  two-dimensional  design  possibilities.  However,  the 
method  is  easy  to  use,  and  good  agreement  with  overall  heat  losses  was  observed 
with  insulated  pipes. 

Numerical  modeling  with  convection  and  radiation  was  demonstrated,  produc¬ 
ing  good  comparisons  of  heat  loss  with  the  experimental  data;  however,  large 
geometries  and/or  large  temperature  differences  across  the  air  gap  were  difficult 
to  model.  The  inclusion  of  radiation  is  required,  and  effects  of  surface  emissivity 
values  can  be  observed.  Numerical  data  were  obtained  only  for  relatively  small 
utilidors,  the  primary  limitation  being  related  to  computational  memory  require¬ 
ments  and  the  matrix  solution  methods.  Future  improvements  in  computational 
methods  and  storage  hardware  may  make  this  analysis  method  practical. 

Numerical  conduction  analysis  using  an  effective  conductivity  produced  rea¬ 
sonable  approximations  to  temperature  distributions  on  inside  surfaces;  the  method 
was  easy  to  use,  and  the  solutions  were  all  obtained  on  a  personal  computer.  Two- 
dimensional  effects  can  be  differentiated,  but  the  information  regarding  tempera¬ 
ture  distribution  within  the  air  gap  is  inaccurate.  The  method  seems  to  be  most 
accurate  for  small  temperature  differences  across  the  air  gap  and  is  relatively 
insensitive  to  minor  changes  in  the  effective  conductivity  value. 


CONCLUSIONS  AND  RECOMMENDATIONS 

Average  heat  losses  can  be  calculated  reliably  for  insulated  pipes  using  one¬ 
dimensional  analysis,  and  these  results  will  compare  well  with  full  (convection 
and  radiation)  numerical  solutions  (of  average  heat  loss). 

A  full  numerical  solution  will  provide  the  best  two-dimensional  analysis.  How¬ 
ever,  the  current  model  may  not  be  able  to  converge  to  a  solution  given  reason¬ 
able  computer  resources.  Ignoring  radiation  in  a  numerical  convection  model  of 
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utilidors  will  have  a  significant  effect  on  the  temperature  distribution  and  will 
result  in  lower  predicted  heat  transfer  rates. 

The  use  of  an  effective  conductivity  for  air  in  a  numerical  conduction  analysis 
will  produce  reasonably  good  temperature  distributions  on  interior  surfaces.  How¬ 
ever,  the  air  temperature  distribution  will  be  in  error.  The  procedure  is  relatively 
insensitive  to  the  effective  conductivity,  the  pipe,  and  enclosure  insulation  domi¬ 
nating  the  heat  loss,  at  least  for  the  cases  investigated  in  this  work. 

A  more  efficient  and  robust  numerical  modeling  approach  is  required.  Two  pos¬ 
sible  improvements  are  to  (1)  convert  the  solution  procedure  to  a  "segregated 
method"  where  each  of  the  governing  equations  are  solved  individually  and  a  Pois¬ 
son  equation  is  substituted  for  the  continuity  equation,  or  to  (2)  incorporate  an 
upwinding  scheme  such  as  the  Petrov-Galerkin  method  into  the  element  quadra¬ 
ture  procedure. 

Some  temperature  data  from  actual  utilidors  are  available;  comparisons  with 
these  data  and  numerical  modeling  of  the  entire  soil  mass  should  be  done.  Addi¬ 
tionally,  comparisons  with  an  effective  conductivity  correlation  in  a  transient 
conduction  model,  including  the  soil  mass,  and  compared  with  field  data,  may 
also  produce  some  interesting  results. 
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_ Average  beat _ 

Left  Side _ Top _ Average  Flux  Conductance 

i  Ra  Nu  Ra  Nu  Ra  (Btu/hr  ft)  (Btu/hr  ft2oF) 


Temperature  CF) _  Average  heat 

Insulation,  pipe  Insulation.  EPS  Bottom  Right  Side  Leftside  Top  Average  Flux  Conduct 


APPENDIX  B:  NUMERICAL  RESULTS 


_ Temperature  (°F) _  Average  heat 

Pipe  Inside  Outside  Bottom  Side  Top  Average  Flux  Conductance 

Filename  Pipe  Insul.  EPS  EPS  Dtemp  7>G2  Nu  Ra  Nu  Ra  Nu  Ra  Nu  Ra  (Btu/hrft)  (Btu/hr  ft2oF) 

sq4gr35r  |  10o|  |  79.13|  65 1  20.87|  89.571  8.24|  1.68E+06I  8.501 1.66E+061  10.181 1.54E+06I  9.05  1.62E+06  26.25  0.343 

sq4gr40r  100  ~76.08|  60|  23.92|  88.04|  8.n|  1.94E+06|  8.3 8 1 1.92E+06|  9.98|  1.79E+06|  8.92|  1.88E+06|  29,60  |  0.337 
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|sq4gr05s  |  100|  96.91 1  95 1  3.09 1  98.461  7.50|2.33E+05|  8.04|2.28E+05|  10.05  2.08E+05  8.56  2.22E+05  3.73  0.329 

sq4grl0s - 100  IQT  ~90  6.19  96.91  7.43  4.73E+05  7.92  4.63E+0S  10.12  4.20E+05  8,49  4.51E+05  7,38  0.325 
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sq4gr30s  ~T00  8  L25  70  18.75|  90.62|  7.18|  1.50E+06|  7.47|  1.48E+06|  9.54|1.35E+06|  8.06|  1.44E+06|  21.04  |  0.306 
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Temperature  (°F)  Average  heat 

Pipe  Inside  Outside  Bottom  Side  Top  Average  Flux  Conductance 

Filename  Pipe  Insul.  EPS  EPS  Dtemp  TAVG2  Nu  Ra  Nu  Ra  Nu  Ra  Nu  Ra  (Btu/hrft)  (Btu/hr  fr°F) 
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Temperature  (°F)  Average  heat 

Pipe  Inside  Outside  Bottom  Side  Top  Average  Flux  Conductance 

Filename  Pipe  Insul.  EPS  EPS  Dtemp  TAIG2  Nu  Ra  Nu  Ra  Nu  Ra  Nu  Ra  (Btu/hrft)  (Btu/hr  ft-°F) 


Tj-rf-rocororOrOro 

0  0*000000 


*4,  __  W  Ln 

-'t'o»222l2 


???????? 

WWWWWWWW 

^OrirnChtS>nN 

SO^ONfiOOOrtVIfS 

rdWW<Nrdinvdod 


oooooooor>|c^i^vo 


WWWWWWWW 

fOVlHtONO\W^ 

rs,t'0^rS*Hn'c 

nosHiNW’t'n^ 


t^OscococNVOcoOs 

intNOOflOhiocnoi 


Sto  >n  >ti  >n  >n  in  <n 

o  o  o  o  o  o  o 

wwotiwwwww 

oo^ai^MooOfn 
r-^ososorooovo 
cd  ~  ~  c4  uS  vd  po 


nHrfMisinOsn 

fnON'cirsio>^Ofn^ 

#d  K  h  K  io  ^  so  ^ 


rfinminmininin 

¥¥?????? 

wwwwwwww 

fNHfSlATfTf'Oin 
o  N  ^  rn  tv  rf,  ni 
TfHoirnTfinihOi 


r-r-sovor-osoco 

-hrf-HOOinrnO 

r^vc*ovoviin»n<n 


00  O  CO  ho  r-  Os  ^ 

oooolr^Ko»nKr|^t 


Nmrt^Hfnoovo 
OOOVO^ffSOOOt^ 
^NTfidodd^rn 


m  r*  oo  on  h  n  rt 

ONOOi^io^in^fn 


^rNtnOin^'OTr 

00«fSNNN*HO 

vbocovdosfsvioo 

Osasoor-soNosnTf 


^  in  x  ^  in  ^ 

oi  Mn  ^  n  o 

co  ro  ro  rO  co  CO 

0  0  0  0*0  0 


Os  00  VO  <N  ON  rz 

4.  o\  ~  ^  00  ^ 

~  CO  d  00  Os  ^ 


tt  o  m  m  <n  <n 

000000 

£  |£  £  £  5  £ 

ro  co  *n  00  co 

tj-  co  ro  -t-  t-~  ro 

Tt  ^  ri  rn  rf  16 


c*  >n  h  o  o  n 

co  ©  r*  os 

00  00  vd 


r)1  in  n  n  n  in 

000000 

w  w  w  w  i  w 

00  ro  Tj-  in  O  CO 

ooh^oo^h 

CO  ^  *  (S  CO  »o 


¥  ¥  ?  ¥  ¥ 

w  w  w  w  w  w 

00  00  co  ro  po  so 

*n  co  ^  v©  ©  v© 

't  ^  ri  rn  o  no 


t  Q  't  ro  n  01 

t**  co  Os  vo  co  O 

NO  NO  NO  so 


s  £  £  £  £  £ 

??¥¥¥? 

m  m  m  m  u  m 

x  h  x  \o  n  h 

«  fj-  in  oe  rn  m 

rf  h  oi  cn  vi  K 


os  0  os  o  <n 

in  CN  00  v©  ro  O 

no  no  mi  »n  »n  vs 


^  r) \rf  I  vo  00  lo 


r-  co  »— 1 

**<  st  no  oo 

K  00  oi  O 

VO  w-i  Tt  rf  CO 


CO  VO  <N  ©  tJ*  vo 
00  [0  Ir-*  r~i  lo  OS 


W-H  4-(  «4-4  W_( 

in  in  vn  in  n  n  n 

h  cl  co  ^  in  io  h 

£  :£  rS  r2  r2  rS  rg 

Vj*  Vr  W  Vt  Vt  Vfr  *t 

o'  cr  0*  or  ct1  cr  a* 

to  vi  Vi  vi  Vi  cn  1/5 


bD  bfl  bfl  bD  bfl  bfl 
0  o  in  m  o  mi 

O  ^  (S  fO  t  ^ 
45  X>  J3  Xi  i  £ 


(N  O  OO 
CO  CO  <M 
0*0  0 


VO  O  t-* 

tS  rn  O 

•*•  n  mi 


sf  n  n 

09° 
$  m  $ 

O  N  NO 
h  o  in 
vd  cvi  co 


a  a  $ 

t-i  jvo  vd 


^  m  in 

000 

£1  w  w 

<n  x  o 
vo  os 
«n  W  p4 


Os  *0  ro 
«  St  n 

odd 


st  rf  in  n 
©  ©  ©  © 

W  lS  i 

m  o  tj-  o 

OO  'Sj-  Tf  — j 

ci  x  h  n 

^  7T  5"  v7 
00  vo  *3-  (N 


IV  ^  m  m 

I?  ?  ?  ¥ 
gsss 

in  v  cs  06 


os  lo  lo  lo 


¥  ¥ 
WWW 

o  I-  V 

os  ~  r*- 


VO  js-s  00 

vd  vd  id 


3  jq  jq 
?  ?  ? 
MUM 
v  in  x 

*3*  <N  Os 


?v  vs  o  >n  in  in  mi 
| 

mmmmmmmm 

©p-s^oor^'^rt— ioso 

„  ONNO-sf^OlOO^Os 
g  fsioOi-iriNnstmi 
*X:r-(NVO(NOsC-NOTf 

o5m»— ioovoco»— 'Ost^- 
a^KvoNoiovdinin 


St  ^  m 

fm 

o  00  o  00 
“  o  in  m 
ffn  oi  i-i 

.fc  vn  —  1 n 
g  ro  o 
•=  vd  vo  in 


n  n  n  m  m 

¥  ¥  ¥  ¥  ¥ 

w  w  w  w  w 

CO  i-4  O  1-s  CO 

ri  n  t  vi  vd 


*n  Im  |sn  Uj*  jrf 


•^osNOinr-nJ^stin 
(jrnxNtninvq\ov} 
gdooHcninr^ONH 
5  t  ro  ro  s-  OOsOs 

n  *— <  s— 1  s—1 


—  N  st  *  h  n  m  ^  On  h-  cn  «n 


jm  |Crj  IW  ,  t  ]  ,  ,  ,  . 

rninmninininininininin 

0’-<<nV7o*— itNcorrtnvcr- 

XXjD  X  XjDXjDjQXjDX 
T?  's?  ^  itm  s?  ^  s?  ^  T?  ^  *5*  s? 
o-o-tri  o'crtro'crircro 

tncnuov-itnoocncntntnincn 


61 


Back  to  Contents 


sq4ib85l - 150  ~92T4  ~7lQ5  ~65  17.09  83.59  4.44  7.67E+05  5.53  7.1 6E+05  9.70  5.70E+05  6.23  6.86E+05  18.58  0.297 

sq4ib951 - 150  ~85j02  "6511  55~19lO  75.47  4.23  9.23E+05  5.31  8.60E-H15  9.46  6.83E-H)5  6.00  8.24E+05  19.77  0.282 

s4ibl05l - 150  "7779  ~56!68  45  21.11  67.23  4.03  U0E-»06|  5.09|  1.02E-H)6j  9.15|8.13E+05|  5.76|9.81E+05|  20.71  |  0.268 


8  C 

a  J2 

■a  1 

U  « 


u  & 

od 

CO 

u 

o 

> 

*1 


£ 


cx 

H 


£ 


o 

CQ 


Pu , 

1 

2 

I' 

I 


P 

Cu 

e 


oo 

vo 

<N 

© 

r- 

on 

oo 

r- 

vo 

*o 

CO 

CN 

IN 

IN 

(N 

(N 

CN 

© 

o 

© 

© 

© 

© 

m 

oo 

VO 

(N 

ON 

r^ 

cn 

Tf 

co 

© 

Tf 

r- 

CO 

vo 

00 

On 

TT 

VO 

vo 

vo 

vo 

vo 

o 

¥ 

¥ 

¥ 

¥ 

¥ 

pq 

W 

W 

pq 

W 

W 

o 

CO 

TC 

© 

VO 

vo 

vo 

vo 

On 

ro 

© 

id 

IN 

CO 

vo 

rT 

cn 

CO 

TT 

NO 

VO 

Tf 

CO 

— * 

ON 

r- 

VO 

CO 

vd 

VO 

VO 

»o 

vd 

vd 

Tf 

VO 

vo 

VO 

VO 

vo 

o 

o 

© 

© 

© 

© 

W 

W 

pq 

W 

pq 

00 

CN 

Tf 

© 

VO 

co 

<N 

CN 

TT 

r- 

TT 

*— 1 

IN 

CO 

»d 

— 

co 

On 

CN 

co 

•o 

VO 

ON 

© 

© 

00 

OO 

00 

ON 

On 

On 

od 

VO 

VO 

vo 

vo 

*o 

o 

© 

© 

© 

¥ 

¥ 

W 

W 

PU 

W 

PQ 

PQ 

Tf 

00 

»o 

r- 

© 

oo 

ro 

VO 

r- 

© 

vo 

co 

»d 

— 1 

CN 

N- 

*d 

rT 

VO 

vo 

© 

r- 

vo 

N- 

00 

vo 

m 

© 

OO 

VO 

voi 

»o 

vo 

VO 

N" 

"=± 

uT 

vo" 

"vd" 

— 

o 

© 

© 

© 

© 

© 

+ 

+ 

+ 

+ 

+ 

+ 

W 

W 

W 

W 

W 

w 

o 

ON 

rr 

NO 

ON 

rq 

VO 

ON 

CO 

ON 

ON 

*o 

nL 

CN 

N* 

vd 

t-' 

o" 

ON 

vo 

N- 

TT 

TJ* 

oo 

TT 

CN 

© 

00 

NO 

TT 

rr 

co 

CO 

rT 

” 

vo" 

to* 

3" 

co 

VO 

vo 

vq 

VO 

CN 

VO 

00 

© 

CN 

Tf 

vo’ 

r- 

VO 

vo 

»o 

T}- 

CO 

-T 

“ 

~ 

rT 

oT 

T?" 

IN 

ro 

CO 

rq 

CN 

CN 

i— 3 

cd 

vo 

rT 

ON 

vo 

«o 

“ 

"vo" 

Td" 

VO 

r- 

vo 

vo 

Tf 

CO 

CN 

“ 

©" 

oT 

oT 

CO 

© 

i> 

ON 

ON 

ON 

On 

oo 

vo 

VO 

rT 

OO 

On 

© 

e'¬ 

VO 

vo 

N- 

CO 

CO 

er 

ro" 

co 

i— t 

ON 

CN 

co 

co 

CN 

o 

VO 

o 

CO 

vd 

ON 

CN 

r- 

r- 

vo 

VO 

rr 

rC 

©” 

©" 

pi 

ipl 

00 

oo 

jl 

B 

00 

:»**» 

vo"1 

VO 

voT 

VO 

»o 

VO 

o 

(N 

ro 

''t 

«o 

r2 

r£ 

jO 

x> 

*cr 

"o' 

"S 

S- 

s- 

on 

OT 

Vi 

Vi 

VI 

VJ 

vo 

© 

ON 

rC 

CN 

ro 

NO 

vo 

© 

NO 

*** 

r- 

r» 

rvl 

rsl 

r— < 

ON 

VO 

vo 

ro 

CN 

O 

oo 

r- 

CN 

CN 

CN 

CN 

CN 

i 

CN 

CN 

CN 

CN 

CN 

CN 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

CO 

vo 

ON 

ON 

r** 

Tf 

ON 

Os 

O 

vo 

VO 

On 

On 

>— < 

© 

NO 

ro 

vo 

CN 

N- 

ON 

00 

O 

l-H 

CN 

CO 

vo 

rd 

ON 

vd 

d5 

CN 

oo 

CN 

rm 

VO 

w** 

»o 

*o 

VO 

VO 

»o 

vo 

VO 

vo 

NO 

vo 

VO 

VC*" 

"3* 

© 

+ 

© 

+ 

? 

? 

? 

? 

© 

+ 

© 

+ 

© 

+ 

© 

+ 

© 

+ 

© 

+ 

© 

+ 

w 

W 

w 

pq 

w 

W 

pq 

pq 

pq 

pq 

pq 

pq 

pq 

pq 

i-H 

© 

© 

CO 

Tf 

00 

VO 

r- 

r~- 

vo 

r-; 

vo 

*"* 

ON 

On 

NO 

© 

N- 

VO 

ro 

vo 

VO 

r- 

v-* 

co 

rd 

*o 

ON 

rd 

ON 

i-H 

CN 

'■T 

»-< 

CN 

cd 

VO 

NO 

t" 

r- 

IN 

o 

CO 

© 

CO 

•o 

CO 

cN 

vo 

CN 

w 

ON 

VO 

(N 

ON 

»o 

© 

vo 

© 

r- 

vd 

VO 

VO 

VO 

'Vt 

Tf 

oo 

vo 

vo 

vd 

VO 

vo 

vn 

vo 

VO 

vo 

vo 

VO 

vo 

TT 

VO 

to 

"3" 

© 

© 

© 

? 

? 

© 

+ 

© 

© 

? 

? 

¥ 

¥ 

¥ 

¥ 

i3 

w 

jq 

8 

W 

tu 

pq 

uq 

pq 

pq 

uq 

pq 

pq 

pq 

oo 

NO 

CO 

r- 

CN 

CO 

00 

ro 

© 

oo 

vo 

CN 

VO 

CN 

oo 

ON 

© 

*o 

ro 

Tf 

tj* 

vo 

vo 

vo 

CN 

rd 

Tj- 

r- 

CN 

ON 

IN 

N* 

— ■ 

IN 

cd 

vo 

NO 

N- 

r- 

oo 

ON 

© 

VO 

Ov 

f-N 

rm 

*o 

CN 

CN 

© 

cN 

VO 

vo 

ro 

o 

CO 

00 

ro 

vo 

r- 

CN 

od 

00 

00 

od 

od 

oo 

ON 

On 

00 

od 

r- 

vo 

vd 

tj- 

»o 

»o 

•o 

VO 

NO 

VO 

»0 

vo 

vo 

vo 

VO 

vo 

VO 

? 

o 

+ 

© 

© 

© 

? 

? 

? 

? 

¥ 

© 

¥ 

© 

© 

PJ 

w 

W 

w 

w 

W 

pq 

pq 

uq 

pq 

pq 

pq 

pq 

VO 

a 

vo 

ro 

ON 

oo 

NO 

N- 

Os 

o 

© 

ON 

vo 

CN 

© 

03 

Ph 

© 

vo 

NO 

■vj- 

c-« 

Tf 

VO 

00 

ro 

Tf 

VO 

CO 

ON 

CN 

* 

r-i 

CN 

<o 

vo 

00 

© 

00 

VO 

nj* 

w 

VO 

N- 

vo 

»o 

VO 

00 

r- 

© 

© 

00 

NO 

Tf 

CN 

ON 

Im 

00 

VO 

r-* 

00 

VO 

c^ 

»o 

VO 

N* 

xr 

CO 

o 

r- 

r- 

r- 

vd 

VO 

vo 

»o 

*o 

1 

TrT 

o" 

“ 

vo" 

‘5 

vd" 

3" 

75" 

vd* 

? 

© 

+ 

? 

? 

© 

+ 

© 

+ 

£ 

B 

© 

+ 

¥ 

? 

9 

¥ 

¥ 

¥ 

¥ 

W 

W 

W 

w 

w 

pq 

pq 

pq 

pq 

pq 

w 

pq 

a 

pq 

CN 

00 

VO 

»o 

r- 

00 

u 

Tt 

© 

IN 

CN 

vo 

vq 

rr 

Tf 

VO 

© 

** 

r- 

vq 

TC 

oq 

vq 

vq 

oo 

OO 

rd 

vd 

3. 

On 

o 

cd 

ON 

IN 

CN 

cd 

rmt 

N- 

©" 

© 

© 

? 

a 

JR 

© 

tn" 

rT 

dr 

rT 

© 

oo 

vq 

vq 

rq 

© 

e 

r- 

cN 

On 

vo 

CO 

cd 

cd 

rd 

rd 

o 

*■5 

vd 

vd 

vd 

*d 

>d 

*d 

a 

’b 

Vi 

vo" 

wT 

oT 

“ 

.5 

o' 

rT 

TT 

dr 

TT 

»d" 

CN 

rq 

rq 

ro 

CN 

1 

ro 

00 

rq 

r- 

CN 

vq 

oo 

© 

vd 

<N 

©* 

CN 

© 

od 

5 

r^ 

od 

© 

CO 

co 

CN 

CN 

i-H 

i 

& 

rj- 

ro 

ON 

On 

On 

oo 

(M 

o 

f-H 

ro" 

rH 

3" 

TT 

rT 

TT 

O 

NO* 

oT 

cd" 

oT 

Td" 

CN 

CO 

rq 

rq 

co 

CN 

CN 

CO 

»o 

»q 

rq 

ON 

<N 

N- 

vd 

© 

> 

cd 

© 

cd 

© 

od 

»d 

1 

*S 

CN 

*— < 

CN 

TT 

o” 

TT 

TT 

©" 

i 

vo 

xdr 

~ 

Id" 

TT 

~ 

TT 

co 

ro 

CN 

Tt" 

PO 

CN 

On 

00 

VO 

i 

u 

* 

t— 1 

*"*  i 

r-^ 

V 

a 

TT 

o" 

oo" 

_ 

VO 

oT 

Cb 

©" 

CO 

N- 

f— < 

[59 

gl 

P3 

E3 

vq 

CN 

CN 

vq 

vq 

CN 

JS 

vq 

VO 

Cl 

KJ 

VO 

vd 

i-H 

CN 

C^ 

od 

vd 

S3 

vo 

On 

CN 

*d 

53 

Ri 

fj 

ro 

co 

CN 

i 

■Vt 

co 

ro 

CN 

00 

id 

’~l 

mm 

HI 

IH 

7T 

TZ" 

” 

rd" 

” 

“ 

r- 

■ 

■ 

\m 

■ 

■ 

■ 

oo 

vq 

*o 

© 

00 

ON 

n 

■ 

m 

■ 

■ 

■ 

vd 

rd 

vd 

rd 

vd 

IN 

■ 

m 

■ 

■ 

■ 

■ 

CO 

co 

<N 

CN 

£ 

1 

■ 

■ 

■ 

©" 

©" 

©" 

©" 

©" 

o’ 

a 

© 

113 

IE3 

© 

© 

© 

m 

N- 

N- 

N- 

N* 

N- 

vo 

VO 

VO 

vo 

© 

© 

© 

© 

B 

© 

e 

vo 

© 

•2 

© 

3 

© 

cN 

r2 

vd 

Is 

vd 

ro 

-2 

© 

*o 

jO 

i 

fH 

X 

ca 

© 

T3 

© 

VO 

3 

© 

VO 

(N 

a 

cd 

»o 

ro 

T3 

-o 

1 

x > 

vo 

3 

X) 

vo 

x> 

»0 

CO 

*o 

Vt 

"s 

'T 

Vt 

dr 

Cn 

CN 

CN 

CN 

CN 

IN 

CN 

IN 

XT 

ar 

XT 

XT 

XT 

l 

cr 

cr 

cr 

cr 

cr 

a- 

cr 

CT 

Vi 

V) 

Vi 

Vi 

Vi 

Vi 

fH 

Vi 

Vi 

Vi 

Vi 

V) 

Vi 

OT 

OT 

VO 

VO 

© 

6© 

VO 

*o 

CN 

T— 4 

© 

© 

© 

© 

rT 

VO 

cn 

»o 

rq 

vo 

CN 

•d 

VO 

vo 

VO 

vo 

© 

+ 

© 

+ 

¥ 

¥ 

pq 

pq 

pq 

pq 

IN 

oo 

On 

VO 

r- 

f— ^ 

vo 

vd 

*-* 

cn 

Tf 

© 

On 

CN 

VO 

ON 

vo 

»d 

*d 

<o 

»o 

VO 

VO 

vo 

¥ 

¥ 

¥ 

¥ 

pq 

pq 

pq 

pq 

IN 

VO 

© 

CO 

r- 

© 

vo 

vd 

— 

cn 

<<t 

*— > 

© 

oo 

VO 

ON 

On 

o 

VO 

IN 

VO 

VO 

vo 

vo 

VO 

© 

VO 

o 

© 

© 

¥ 

+ 

+ 

+ 

pq 

pq 

pq 

pq 

CN 

1—4 

r- 

m 

r-* 

00 

— H 

vo 

vd 

-H 

cn 

Tf 

© 

Tf 

vo 

*** 

CN 

VO 

CN 

© 

VO 

*o 

«o 

vd 

o" 

” 

“ 

To" 

¥ 

¥ 

¥ 

9 

a 

a 

w 

ON 

pq 

vo 

oo 

00 

vq 

vd 

1—1 

cd 

tT 

ON 

— M 

cn 

vo 

•-H 

ON 

»d 

vd 

vd 

Tf 

dT 

"•d" 

ON* 

CN 

N; 

vo 

rq 

oo 

Tf 

© 

vd 

r- 

“ 

r- 

VO 

rT 

o’ 

rn" 

vq 

© 

tq 

cd 

— 3 

od 

vo* 

«— i 

mrn 

CN 

~ 

tn 

vd" 

vd" 

r- 

vo 

VO 

N* 

»— < 

oo" 

©" 

r- 

tT 

ON 

rq 

vq 

vd 

od 

rm 

cd 

r- 

vo 

VO 

*o 

1 

1 

1 

1 

© 

© 

© 

IK=9 

oo 

00 

00 

oo 

o 

o 

o 

o 

vo 

vo 

vo 

vo 

© 

a 

a 

§ 

IN 

3 

IN 

cr 

cr 

cr 

cr 

OT 

OT 

OT 

OT 

cn 

ON 

r- 

vo 

cn 

CN 

IN 

V* 

© 

o 

© 

© 

CN 

© 

ON 

CN 

00 

F“! 

CN 

vo 

cN 

*o 

vo 

vo 

¥ 

f 

¥ 

¥ 

pq 

pq 

pq 

pq 

ON 

N* 

rr 

© 

rr 

r- 

00 

oo 

IN 

Tf* 

r- 

oo 

N- 

© 

© 

vo 

oo 

VO 

cn 

vd 

tT 

*o 

Vp 

vo 

vo 

¥ 

9 

¥ 

¥ 

pq 

pq 

pq 

pq 

g 

VO 

VO 

cn 

r* 

£ 

oo 

CN 

VO 

cn 

oo 

vo 

r- 

© 

»o 

On 

VO 

rT 

VO 

TT 

vo 

VO 

VO 

vo 

¥ 

¥ 

¥ 

¥ 

pq 

pq 

pq 

pq 

vo 

oo 

00 

r-- 

vo 

r» 

oo 

od 

CN 

'•t 

r^ 

00 

CN 

© 

CN 

CN 

vo 

m 

*o 

Tf 

^r 

vo" 

¥ 

© 

+ 

¥ 

9 

a 

pq 

w 

pq 

© 

© 

00 

oo 

oo 

ON 

rm 

od 

CN 

tT 

rT 

on" 

00 

VO 

N- 

a 

© 

nf* 

N" 

__T 

"vo* 

rT" 

©* 

CN 

CN 

rq 

oo 

© 

vd 

cn 

cn 

cn 

CN 

-H 

rm 

to 

rT 

rq 

vq 

nt 

cd 

mm 

On 

rT 

mm 

—i 

CN 

vn" 

Id" 

~ 

*d" 

cn 

CN 

_ i 

On" 

ON 

oo" 

CN 

vo 

vq 

vd 

od 

© 

CN 

cn 

CN 

IN 

1 

II 

II 

1 

© 

nr 

© 

9 

1 

•o 

■o 

•T3 

T3 

vo 

»o 

vo 

vo 

g 

5 

!3 

cn 

*CJ 

IN 

IN 

IN 

CN 

cr 

cr 

or 

cr 

OT 

OT 

OT 

OT 

Go  to  Contents 


62 


_ Temperature  (°F) _  _ Average  heat 

Pipe  Inside  Outside  Bottom  Side  Top  Average  Flux  Conductance 

Filename  Pipe  Insul.  EPS  EPS  Ptemp  TAPV2  Nu  Ra  Nu  Ra  Nu  Ra  Nu  Ra  (Btu/hrft)  (Btu/hrft2oF) 
1-  ft  x  1-ft  enclosure,  2.375-inch  pipe,  emissivity  of  pipe:  0.9,  cmissivity  of  EPS:  0.6 _ 
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Temperature  (°F)  Average  heat 

Pipe  Inside  Outside  Bottom  Side  Top  Average  Flux  Conductance 

Filename  Pipe  Insul.  EPS  EPS  Ptemp  T AVg  Nu  Ra  Ku  Ra  Nu  Ra  Nu  Ra  (Btu/hrft)JBtu/hrft^ 
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Temperature  (°F)  _ Average  heat 

Pipe  Inside  Outside  Bottom  Side  Top  Average  Flux  Conductance 

Filename  Pipe  Insul.  EPS  EPS  Dtemp  7>c2  Nu  Ra  Nu  Ra  Nu  Ra  Nu  Ra  (Btu/hrft)  (Btu/hr  ft2oF) 
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sq2ic051  |  200ll96.48|  19S.65|  195|  0.83|  196.07|  10.301 3. 22E+04I  11.171 3.1 1E+04I  13.23I2.88E+04  11.62  3.06E+04  1.98  0.648 

sq2icl51  200  189.29  186.87  185  2.41  188.08  10.00  9.87E+04  10.93  9.52E+04  13.43  8.67E+04  11.43  9.33E+04  5.59  0.632 

sq2ic251 - 20ol  182.01 1778.03|  175|'  3.98|  180.02|  9.7o|l.72E+05|  10.60 [  1.66E+05|  13.33|  i.50E+05|  11.14|  1.62E+05|  8.89  |  0.609 


_ Temperature  (°F) _  Average  heat 

Pipe  Inside  Outside  Bottom  Side  Top  Average  Flux  Conductance 

Filename  Pipe  Insul.  EPS  EPS  Ptemp  T AVg  Nu  Ra  Nu  Ra  Nu  Ra  Nu  Ra  (Btu/hrft)  (Btu/hr  ft;°F) 

i-ft  x  1-ft  enclosure,  2.375-inch  pipe  with  1  inch  of  insulation,  emissivity  of  pipe  insulation:  0-5,  emissivity  of  EPS:  0.9  _ _ _ 
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